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We classify all possible 36 gap-opening instabilities in graphene-like structures in two dimensions, 
i.e., masses of Dirac Hamiltonian when the spin, valley, and superconducting channels are included. 
These 36 order parameters break up into 56 possible quintuplets of masses that add in quadrature, 
and hence do not compete and thus can coexist. There is additionally a 6th competing mass, the one 
added by Haldane to obtain the quantum Hall effect in graphene without magnetic fields, that breaks 
time-reversal symmetry and competes with all other masses in any of the quintuplets. Topological 
defects in these 5-dimensional order parameters can generically bind excitations with fractionalized 
quantum numbers. The problem simplifies greatly if we consider spin-rotation invariant systems 
without superconductivity. In such simplified systems, the possible masses are only 4 and corre- 
spond to the Kekule dimerization pattern, the staggered chemical potential, and the Haldane mass. 
Vortices in the Kekule pattern are topological defects that have Abelian fractional statistics in the 
presence of the Haldane term. We calculate the statistical angle by integrating out the massive 
fermions and constructing the effective field theory for the system. Finally, we discuss how one 
can have generically non-Landau-Ginzburg-type transitions, with direct transitions between phases 
characterized by distinct order parameters. 



I. INTRODUCTION 

Many of the physical properties of graphene are cap- 
tured by a one-band tight-binding electronic Hamiltonian 
with uniform, real-valued, and nearest-neighbor hopping 
amplitude whereby: (i) electron-electron interactions arc 
ignored; (ii) spin-orbit interactions are ignored; (iii) the 
electronic band structure is replaced by two conical dis- 
persions centered about two non-equivalent points, the 
Dirac points, in the first Brillouin zone; and (iv) the 
coupling to electro-magnetic external fields is governed 
by the minimal substitution. For instance, graphene dis- 
plays an integer quantum HaJieffcct (IQHE) as a function 
of the applied bias Y|)ltage,Eia and it shows an universal 
optical conductivityjj Both these properties can be un- 
derstood within the non-interacting electron picture. 

Although most experiments observe the massless Dirac 
spectrum assumed in (iii), electronic instabilities in the 
form of single-particle spectral gaps (mass gaps in short) 
can be triggered by externp|l perturbations such as some 
commensurate substrates,cl or large enough magnetic 
fields that can change thei-balance between the kinetic 
and the potential energy.BQlD In this paper we study a 
number of issues pertaining to Dirac fermions in two- 
dimensions when a mass gap is opened in the fermionic 
spectrum by different non- vanishing order parameters. In 
particular, we shall study in great detail the simpler case 
when there is no superconducting instabilities and spin- 
rotation invariance is maintained, in which case there are 
only 4 possible masses. We derive in this simpler case the 
effective action when the massive fermions are integrated 
out, and read from this action the fractional statistics of 
topological defects in the mass order parameters. We also 
present a complete classification of all possible masses (36 



in total) in the general case where any spin, valley, and 
superconducting instabilities are permitted. 

In the simpler spinless problem (or, more realistically, 
the problem when spin-rotation invariance is never bro- 
ken), the 4 different masses that can be added to the two- 
dimensional Dirac equation representing graphene are 
the following. One perturbation is a staggered chemical 
potential, taking values +(1^ and —fj.^ in the two sublat- 
tices of the honeycomb lattice of graphene say. It opens 
a gap 2|^g| at the two Dirac points.! A second mass gap 
2 1 77 1 arises by adding directed next-nearest-neighbor hop- 
ping amplitudes in the presence of fluxes, but such that 
no net magnetic flux threads a hexagonal Wigner-Seitz 
unit cell of graphene say. XJris perturbation breaks time- 
reversal symmetry (TRS).a Finally, a real- valued mod- 
ulation of the nearest-neighbor hopping amplitude with 
a wave vector connecting the two Dirac points (i.e., a 
Kekule diaaerization pattern for graphene) also opens a 
gap 2|A|.EiJ This real- valued modulation of the nearest- 
neighbor hoppings is parametrized by the complex or- 
der parameter A = Re A -I- ilm A whose phase controls 
the angles of the dimerization pattern. This mass cor- 
responds to two real masses Re A and ImA, bringing 
the total number of real-valued masses that conserve the 
electron number and spin-rotation symmetry (SRS) to 
four. 

If the order parameters 77, and A are not uniform, 
but vary in space and contain topological textures, then 
midgap states in the massive Dirac spectrum can appear. 
Exarnpjes are static line defects at which and 77 change 
signSjEil and static point defects represented by vortices in 
the phase of A.I12I As occurS[ai;|a,static domain wall in one- 
dimensional polyacetylene,ll3E3L3 a fractional electronic 
charge is exponentially localized in the vicinity of a static 
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charge ±1 vortex in the phase of A.0 

The value of the fractional charge that is bound to a 
vortex in the phase of A also depends on whether the vor- 
tex is drcsacd_with a half flux of the axial vector potential 
or not.Ej'EEl When the axial gauge flux is absent (log- 



arithmically confined case 
be tuned continuously as 
where m -^/jAp + fi^. 



the value of the charge can 
Junction of the ratio fJ-s/m 
It is independent of the 
ratio A's/"^ when the axial gauge half flux is present (de- 
conflncd case) , for thcj_diarge is then pinned to the ratio- 
nal values Q — ±l/2.E3ll3 These values of the fractional 
charges persist as long as the magnitude of the 
breaking mass \ri\ is smaller than the mass scale m 
There is a phase transition at |?7| = m. For I77I >pHi_the 
fractional charge bound to the vortices vanishes. EJM 

Just like the charge, the statistical phase Q acquired 
upon the exchange of two vortices depends on whether 
the vortex in the phase of A is screened or not by the 
axial gauge flux. In this paper, we derive the statistical 
angle from the effective action obtained upon integrating 
out the massive fcrmions. (We thcmbv-ctisolvc conflicting 
claims about in the literature.lijO'E£l) The statistical 
angle depends on the interplay between the magnitude 
of the TRS-breaking mass rj and the magnitude m of the 
TRS masses. There arc phase transitions at the lines 
|?7| = TO depicted in Fig. |l| that separates regions domi- 
nated by the TRS-breaking masses and those dominated 
by the TRS-preserving mass 77. The statistics @ jumps for 
both the screened and unscreened vortices at the phase 
boundaries. 

When unit vortices in A_are screened by an axial gauge 
flux, they are deconflncdo Their statistics is well-defined 
in a dynamical sense and it takes universal values inde- 
pendent of the ratio /.(^ on both sides of the transition. 
We show that 



G = when m > \r]\ 



and that 



— = sgn(77) ( 

TT 



sgnyy 



when I77I > TO. 



(1.1a) 



(1.1b) 



Along the lines = m in the zero-temperature phase 
diagram of Fig. ^, the gap in the Dirac spectrum vanishes. 
At criticality, the notion of point-particles is moot and 
so is the question of their quantum numbers. 

A remarkable complementarity has emerged. Defects 
carry either a fractional charge Q = ±1/2 but no frac- 
tional statistical phase when the breaking of TRS is not 
too strong {\ri\ < m), or no fractional charge but a frac- 
tional statistics O/tt = ±1/4 when the breaking of TRS 
is dominant {\ri\ > to). 

When unit vortices in the order parameter A are not 
accompanied by ^ axial gauge flux, they are logarith- 
mically confined.EHl Although their statistics is not well- 
defined dynamically, it is nevertheless possible to create 
them and exchange them by external means. If so, both 
their charges and statistics acquire a dependence on all 



masses 77, fi^, and A, that we compute analytically and 
test numerically here in this paper. 

We then go beyond the simpler spinless case with only 
4 masses, and we classify all 36 masses in the general case 
where any spin, valley, and superconducting instabilities 
arc allowed. These 36 order parameters break up into 
56 possible quintuplets of masses that add in quadra- 
ture (to a value to^), and thus do not compete with one 
another. The Haldanc mass, the generalization of the 
T] mass above, competes with all the other 35 masses, 
and thus one has generically a quantum phase transition 
when 1 77 1 = TO. We argue that these 5-tuplets provide a 
rich playground for Landau-forbidden continuous phase 
transitions. Wc discuss in the paper how any U(l) or- 
der parameter in a 5-tuplet can be assigned a conserved 
charge and supports topological defects in the form of 
vortices. A pair of U(l) order parameters in a 5-tuplet 
is said to be dual if the vortices of one order parame- 
ter binds the charge of the other order parameter and 
vice versa. A continuous phase transition can then con- 
nect directly the two dual U(l) ordered phases through 
a confining-dcconflning transition of their vortices. 

This paper is organized as follows. We deflne the 
relevant continuum Dirac Hamiltonian and review its 
symmetries for the simpler problem with only 4 masses, 
that encodes the competition between a charge-density, 
a bond-density, and an integer-quantum-Hall instability 
at the Dirac (charge neutral) point of any graphene- 
like two-dimensional electronic system in Sec. We re- 
veal a hidden non-Abelian structure of the field theory 
in Sec. Ill that plays an important role when deriving 
the charge and statistics of quasiparticles. The fermions 
are integrated in the background of these 4 order pa- 
rameters and of the U(l)xU(l) gauge fields to leading 
order in a gradient expansion in Sec. IV. The effective 
low-energy and long-wave length interacting field theory 
thereby obtained is a Anderson-Higgs-Chern-Simons field 
theory for bosonic fields: two U(l) gauge fields and one 
phase field. The induced fractional fermion number and 
the induced fractional Abclian statistical phase in the 
Anderson-Higgs-Chcrn-Simons field theory of Sec. IV are 
computed in Sec. [v| and Sec. VI, respectively. The nu- 
merical calculation of the fractional charges and statisti- 
cal phases within a single-particle (mean-field) approxi- 
mation that viol ates the U(l)xU(l) gauge symmetry is 
presented in Sec. VII. A microscopic (lattice) model shar- 
ing the same U(l)xU(l) gauge symmetry and low-energy 
long- wave-length particle content as the Anders on-Hi ggs- 
Chern-Simons field theory is constructed in Sec. VIII . Ei- 
ther by en largin g the particle content of the lattice model 
from Sec. VIII or by allowing additional magnetic, spin- 
orbit, or superconductivity instabilities to compete with 
the charge-density, bond-density, and integer-quantum- 
Hall instabilities in graphene-like two-dimensional sys- 



tems, we are lead to a classification presented in Sec. IX 
of all 36 competing orders of a Dirac Hamiltonian rep- 
resented by 16-dimensional Dirac matrices that encodes 
the quantum dynamics of electrons constrained to a two- 
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FIG. 1: Phase diagram parametrized by the TRS mass m and 
the TRS-breaking mass 77. There are three regions delimited 
by the boundaries = m, in each of which the spectral 
gap does not close. The boundaries \t]\ — m are lines of 
critical points at which the spectral gap closes. When the 
vortices are screened by half of an axial gauge flux, they carry 
a fractional fermionic charge of \Q\ = 1/2 with the vanishing 
statistical phase = under pairwise exchange in regions 
for which TRS is weakly broken, i.e., the painted region m > 
I77I. Unit vortices are charge neutral but acquire the non- 
vanishing statistical phase |0| = 7r/4 under pairwise exchange 
in reg ions for which TRS is strongly broken, i.e., \ri\ > m. [See 
Eqs. (1.1).] When the vortices are not screened by the axial 
gauge flux, the charge Q acquires a dependence on the ratio of 
the chemical potential and m, for \r]\ < m, and Q vanishes 
for \rj\ > m. The statistics also depend on which phase one 
sits, but it is non-zero for any 77 7^ 0, and it i s rel ated to the 
value of the charge, as shown in Sees. VI and VU. 



dimensional space, as occurs in graphcnc at the charge 
neutral point say. We conclude in Sec. ^ and relegate 
some intermediary steps to the Appendix. 



II. HAMILTONIAN AND SYMMETRIES: 
SPINLESS CASE WITH 4 MASSES 

The continuum model under consideration in this pa- 
per is defined by the sccpad-quantizcd planar Hamilto- 
nian H := J d'^rH whcrcEl 



n 



U 



scalar 7 



Ti-Q + ''^gaugc 

- i/j'^a ■ (—id) ijj, 
-gauge := V''^" • (0 + 0575)^, 

The 4 components of the spinor- valued operator 



(2.1a) 



(2.1b) 



obey the equal-time fcrmion algebra 
{^,(r),^;(r')} = 



(2.1c) 



The representation (2.1b) is here fixed by the indices A 
and B that distinguish the two triangular sublatticcs of 
the honeycomb lattice and the indices -I- and — that 
distinguish the two inequivalent Dirac points (valleys) 
of graphene. With this choice, the 4 Dirac matrices 
= a^, = o? , ot = c? = i?, and /3— arc defined 
by their 4-dimensional chiral reprcscntationE2l 



a 



T 

-T 

-T, 



0-3 ®r.^ = R, 



(2.1d) 



'0' 



where the 2x2 unit matrix Tq and the three Pauli ma- 
trices , Tj, and T3 act on the sublattices indices (A, B) 
while the 2x2 unit matrix ctq and the three Pauli matri- 
ces (Ti, cr2i (73 act on the valley indices (-|-, — ). The 
matrix 



75 =7^ 







(2.1e) 



acts trivially on the sublattices indices while it acts non- 
trivially on the valley indices, i.e., (l±75)/2 is a projector 
on the + and — valley indices, respectively. In (3+1)- 
dimensional space and time quantum electrodynamics, 
the eigenspaces of (1 ± 75)72 define the chiral indices, a 
terminology that we shall also use in this paper. The 
external (background) real- valued fields a = (ai,a2), 
— (ctsi, 052)7 Q = — argA, /ig, and r\ are space- 
and time-dependent fields. Their microscopic interpreta- 
tion is the following. 

A strong uniform magnetic field (rptational of a) is re- 
sponsible for the IQHE in grapheneO A vector field 
leaeades changes in the curvature (ripples) of graphene, 
E3c3 and it can also encojdfi-defective coordination num- 
bers at apical defects. t3oE3 A constant /ig realizes 
in graphene a staggered chemical potential and opens 
an electronic spectral gap .13 A constant 77 realizes in 
graphene a directed next-nearest-neighbor hopping am- 
plitude without net magnetic flux through the Wigner- 
Seitz cell of the honeyconib lattice and it also opens 
an electronic spectral gap.B A constant A realizes in 
graphene a Kekule distortion of the nearest-neighbor 
hopping aipplitude and, again, opens an electronic spec- 
tral gap£3 The 4 space- and timc-indcpendcnt Re A, 
ImA, and 77 exhaust all possible ways for the open- 
ing of a spectral gap in the single-particle spectrum of the 
kinetic Dirac kernel a. ■ {—id), as /3, /?75, R, and ia^a^ 
generate the largest set of traceless and Hermitian 4x4 
matrices that anticommutes with a- (—id). The 3 masses 
Re A, ImA, and jig are compatible, i.e., they open the 
gap 2m, where 



(2.2) 



for /?, Pj^, and R anticommute pairwise. On the other 
hand, the mass rj competes with the mass m, as ia]a2 
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commutes with (3, P"f^, and R (the competition between 
r] and m leads to a phase transition when I77I = m, which 
shall be important in the discussion of fractional statistics 
in this paper). The fields a, a^, A, and 77 have also ap- 
peared in the context of (a) slave-boson treatments of an- 
tiferromagnetic spin-1/2 HdpiJsjipg model on the square 
lattice in the tt-Aux phase j23Ej'L2lEil and (b) Anderson lo- 
calization for electrons hopping on a square lattice with 
a flux of half a magnetic flux quantum-iprji-plaquette, i.e, 
the square lattice with 7r-flux phase.ooEj 



A. Symmetries 



The model defined in Eq. possesses a number of 
symmetry operations that we list below and utilize in the 
paper. 



1. Time-reversal symmetry 

In the Heisenberg representation, 

H{t) H{~t) (2.3) 
under the anti-unitary transformation 

4,{r,t) ^ (tK^,^ {r,-t), 

a(r, t) -a{r, -t), T]{r, t) -r]{r, -t), ^2A) 
a^{r, t) ^ a^{r, -t), 0{r, t) ^ 9{r, ~t), 
|A|(r,i)^ |A|(r,-t), fi^{r , t) ^ fi,{r , -t) , 

where complex conjugation is represented by K and 



T Pa'-,, 



T, 

ri 



(2.5) 



is a unitary, Hermiti an (and thus symmetric) ma- 
trix. Transformation (2.4) realizes reversal of time in 
graphene, for T exchanges the two valleys while acting 
trivially on the sublattice indices. Moreover, transfor- 
mation ( ^.4[ ) realizes reversal of time for an effectively 
spinless single particle, for T is symmetric. Hamiltonian 
H is time-reversal symmetric^and can be represented by 
real- valued matrix elementSjEj if all background fields are 
static while 



?7 = 0. 



(2.6) 



2. Sublattice symmetry 

Always in the Heisenberg representation, 

H{t) ^ ~H{t) (2.7) 



under the unitary transformation 

a,(r,t) - 
|A|(r,t) 



a{r,t), 7j{r,t) ^ -7]{r,t), 
-.|A|(r,i), fi^{r,t)^ -fi^{r,t), 



(2.8) 



where 



R :-- 



- I '3 









CT3 (g) T3 = i?* 



(2.9) 



is a diagonal, unitary, and Hermitian matrix. Transfor- 
mation {2.i) realizes in graphene the change of sign of 



the single-particle wave functions on every sites of the 
honeycomb lattice belonging to one and only one trian- 
gular sublattice. The single-particle eigenstates of the 
conserved Hamiltonian H obey the spectral symmetry 
(SLS) by which any single-particle eigenstate |^') with a 
non-vanishing energy eigenvalue e has the mirror eigen- 
state Rl'if) with the non- vanishing energy eigenvalue — e, 
if all background fields are static while 



Ms = ^7 = 0. 



3. Continuous gauge symmetries 



(2.10) 



We now turn to the continuous symmetries obeyed by 
the Dirac Hamiltonian (pTl) in the Heisenberg represen- 
tation. To this end, we make use of 



0=[75,a] = {75,/3} = [75,fi]- 



(2.11) 



The commutators and anticommutator (2.11) imply that 



Hit) ^ Hit) 
under the U(1)(8)U(1) local gauge transformation 



lAI 



(0+^575)^^ 
^5 - 



a 



Ms ^ Ms 



^ a d(j), 



(2.12) 



(2.13) 



generated by the two space- and time-dependent real- 
valued smooth functions (j) and (f>,. The microscopic ori- 
gin of the global U(l) gauge symmetry generated by 4> 
is conservation of the electron number in graphene. For 
planar graphene, the continuous global axial U(l) gauge 
symmetry generated by 0g is broken as soon as the cur- 
vature of the tight-binding dispersion is accounted for 
so that the Dirac points are not anymore decoupled. We 
shall nevertheless impose the local axial U(l) gauge sym- 
metry at the level of the approximation captured by the 
Dirac Hamiltonian ( 2A ) and see through its consequences 
in this paper. (We do provide a microscopic example of 
a lattice model that realizes the local axial U(l) gauge 
symmetry in Sec. VHl.) 
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III. PATH INTEGRAL FORMULATION OF 
THE MODEL WITH 4 MASSES 

For our purposes, it will be more convenient to trade 
the operator formalism for an effective partition func- 
tion defined by integrating over the Dirac fermions in 
the background of the gauge fields a and and of the 
scalar fields A, fi^, and rj. We will demand that this ef- 
fective theor y captures the U(1)(8)U(1) local gauge sym- 
metry ( |2.^ ). This is possible in odd-dimensional space 
and time,E3 for the Grassmann measure can be regular- 
ized without breaking the U(1)(8)U(1) local gauge sym- 
metry of the Lagrangian. Of course, maintaining the 
U(1)(8)U(1) local gauge symmetry can only be achieved 
if the phase 9 — — arg A of the Kekulc background field 
A is also included as a dynamical field. For simplic- 
ity but without loss of generality as far as the computa- 
tion of the charge quantum number and statistical phase 
are concerned, the masses m and rj will be taken to be 
space- and time-independent parameters, while A and /ig 
vary in space and time (with m = y^Ap+T'f constant) 
through 9 = —arg A and cos a = Hg/m. (For simplicity, 
we shall also focus on the case where fj.^ is also constant 
in space and time, with the exception of near the vortex 
core, where A — > 0, so fj.^ has to adjust as to keep m 
constant.) 

Thus, we seek the effective field theory defined by the 
Grassmann path integral 



C. 



(3.1a) 



where we have also included the time-components Cq 
(TRS but SLS breaking) and aso (SLS but TRS break- 
ing) of the U(1)(8)U(1) gauge fields to maintain space and 
time covariance. The independent Grassmann-valued 
fields over which the path integral is performed are the 
4-components spinors ip and t/i. They depend on the 
contravariant 3-vectors x^^ ~ {tii') [covariant 3-vectors 
= (t, — r)] and we will use the repeated summation 
convention x^y^ = x'^ip — x^y^ — x^y^. We have defined 
the four gamma matrices 

7° := /3, 7I := /3a\ 7^ := fSa"^ 7^ := /3a^ (3.1b) 

for which lowering and raising of the greek indices 
/i, = 0, 1, 2 is achieved with the Lorentz metric = 
diag(l, — 1, — 1). The 4 matrices 7", 7^, 7^, and 7"^ obey 
the usual Clifford algebra in Minkowsky space in the chi- 
ral representation, i.e., 75 = i7''7^7^7'^ is diagonal. We 
have also defined the matrix 



Ml 



m,r, ■= m ("iMi + 712^/2 + iT-sM^) + '77.57 



1, 



Mo 



-«75: 



/3a 



(3.1c) 



for which we do not distinguish upper and lower latin 
indices a, b — 1,2,3 as they are contracted with the 



Euclidean metric S^^ = diag(l,l,l). (Notice that be- 
cause space and time is (2-1-1) dimensional, we can use 
the gamma matrix 7^^ to open a spectral gap by taking 
M3 = 7^^.) The space and time dependencies in M„ ,^ fol- 
low entirely from those of the phase arg A. Indeed, while 
the masses 77 and m are constant in space and time, the 
direction of the unit vector n with the 3 components 



I A| cos 61 



I A| sin 61 



^ (3.1d) 

TO 



can vary in space and time. 

The U(1)(8)U(1) local gauge symmetry (|2.13|) h as be- 
come the invariance of the Lagrangian in Eg. ( |3.l| ) under 
the U(1)(8)U(1) local gauge transformation 



a 



2'/'5- 



d. 



></^5> 



(3.2) 



In spite of appearances [ipip ■(/'exp(2ic/)575) t/i], the 
Grassmann Jacobian i nduc ed by the U(l)(g)U(l) local 
gauge transformation ( ^.13 ) is unity and does not pro- 
duce a quantum anomaly in (2-1-1) dimensions (odd 
space-time dimension) £3 

We take advantage of the fact that and -0 are in- 
dependent Grassmann integration variables to bring the 
algebra obeyed by the 6 matrices 7^ = 0, 1, 2 and 
a = 1, 2, 3 to a form that will simplify greatly the eval- 
uation of the partition function (3.1). Under the non- 
unitary change of integration variable 



V' =: X757^ 



V' -■ X, 



(3.3) 



the partition function (pj) becomes 



(3.4a) 



-C™,, = X (r^*a^ + r^S^, - mn^E, - r,) X, 
where the matrices 

r^:=757V, Sa:= 757'Ma, (3.4b) 

obey 

{^^ F'^} = 2.9'-, [S,, E,,] = ze,,,I]„ [F", EJ = 0, 

(3.4c) 

for fi,!/ = 0,1,2 and a, b,c = 1,2,3 and we have re- 
grouped the gauge fields into 



(3.4d) 



following the prescription 



bl:=hl:=Q, bl :^ +a,^. (3.4e) 



Notice that 



S3 = -75 



(3.5) 
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so that the sym metr y under the U(1)®U(1) local gauge 
transformation (3.2) has become the invariance of the 
Lagrangian in Eq. (3.4) under 



h]. 



(3.6) 



A. Hidden U(2) non-Abelian structure 

To make the U(2) non-Abelian structure explicit, ob- 
serve first that the mass mn^Y,^ is an element of an su(2) 
Lie algebra. Indeed, there exists a 4 x 4 matrix U repre- 
senting an element of SU(2) generated by Y^^ a = 1,2,3 
such that 



TO n„T., 



mUY^U''. 



(3.7) 



We then infe r that the partition functions (3.1) or, equiv- 
alcntly. (^^) are special cases of the more general parti- 
tion function 



Z := / V[x,x]cxp ( i I (fx£ 



771 , t; 



(3.8a) 



where 

B^{x) = bl{x)+bl{x)j:-, M = 0,l,2, (3.8b) 

are arbitrary elements of the Lie algebra u(2) = u(l) 
su(2) and 

U{x) = e™«(^)e™-(^)^s u„{x),u^{x) G M, (3.8c) 

is an arbitrary element of U(2). As the mapping be- 
tween the unit vector n{x) and U{x) is one to m any , the 
Lagrangian and the Grassmann measure in Eq. (3.8) are 
both invariant under the local U(2) gauge transformation 



B, 



X Vx, 



(3.9a) 



[/-> 

parametrized by the smooth space- and time-dependent 



V{x) 



e U(2), 



(3.9b) 



and under the global U(l)xU(l) transformation 
parametrized by the re al- v alued numbers (/)q and (j)^. 



The transformation (3^) or, e quiva lently, (3^) is rep- 
resent ed b y the tran sformation (3.9) with B^ given in 
Eqs. (3.4cl) and ( ^.4c ) and U given by 



U = e+ies3/2g-iQS2/2g-jes3/2 



(3.10) 



whereby the unit vector (3. Id) is parametrized by 



n ~ (sin a cos 0, — sin a sin 9, cos a) 



(3.11) 



(Recall that cos a Ms/'^i sin a := |A|/to, and that the 
phase 6 = — argA is space and time dependent.) 

A gradient expansion for the partition function (|3^ 
with an arbitrary space and time dependent U G SU(2) 
but with B^^ = and 77 = was performed by Jaroszewicz 
and shown to produce the effective action for the 0(3) 
non-lmftarpSpHWi model (NLSM) modified by a Hopf 
term.L3'L3'L3c2lcil This Hopf term was shown by Chen 
and Wilczek to vanish as soon as the TRS-breaking mass 
77 is larger in magnitude than the TRS mass m. Chen 
and Wilczek also showed that an Abelian Chern-Simons 
term for a non-vanishing 6jJ = is present if and only 
if I77I > m. 

Hopf or Chern-Simons terms can cause the fraction- 
alization of quantum numbers. Although charge frac- 
tionalization can here also be deduced from the presence 
of mi dga p single-particle stateSpXj£_l|hij_Qixac Hamilto- 
nian (2J^) in static backgrounds ,tHlt^LErL3EZl it is natural 
to explore the emergence of fractional statistics under 
the exchange of point-like quasiparticles by exploring the 
fully dynamical theory encoded by the partition func- 
tion (3.4). To this end, it is essential to preserve all sym- 
metries as we did up to now. The point-like quasiparticle 
whose braiding statistics we shall derive are vorticesll2l 
in the dynamical phase 9 = —argA, including the case 
when they are accompanied by axial gauge half flujpes in 
flg^ that screen the interactions between vortices£3 



IV. DERIVATIVE EXPANSION AND THE 
EFFECTIVE ACTION 

It is known that the Dirac Hamiltoniaur 
static backgrounds can support zero modes 
This can be of a nuisance when computing a fermion de- 
terminant. However, it is possible to elegantly dispose of 
this difficulty with the help of the observation made by 
Jaroszewicz that a non-singular U(2) gauge transforma- 
tion on the Dirac Kernel in the partition function (3.8) 
can turn a single-particle midgap state into a single- 



particle thKCshold state without changing the spectral 
asymmetryOo This is achieved by redefining the G rass - 



mann integration variables in the partition function (3.8 
according to 



X =: X 



X 



Ux'- 



The partition function (3.8) becomes 



Zln,n[B'^,] ■■= J X>[x',x']exp (^i J (fxC[„ ,^ , 
X' {r'-id^ + r^S; - mY, - -q) x', 

where 

B' = U^BJJ + UhdJJ 

H' H' M 



(4.1) 



(4.2a) 



(4.2b) 
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need not be a pure gau ge because of the term W B^U . 

The symmetries ( 3.9) of the Lagrangian and the Grass- 
mann measure in Eq. (3.8) become the invarian ce o f the 
Lagrangian and the Grassmann measure in Eq. (4.2) un- 
der the local U(2) gauge symmetry 



X 



X 
B 



(4.3a) 



parametrized by the space- and time-dependent V{x) £ 
U(2) and imdcr the global U(l)xU(l) gauge symmetry 



X 



xw, 



X 



U^UW, (4.3b) 



parametrized by the space and time independent W 
exp(i(/)Q) exp(i03E3). Notice that 

W'^ B 



W 



(4.4) 



under the transformation ( [4.3[ ). 

Evidently, the transformed Dirac fcrmions arc local 
U(2) gauge singlets. Thus, by dressing the original 
Dirac fermions into local U(2) gauge singlets, any midgap 
single-particle states from the original static Dirac Hamil- 
tonian has migrated to the threshold of the continuum 
part of the transformed single-particle spectrum, pro- 
vided the single-particle spectral gap has not closed, i.e., 
m ^ I77I in the parameter space {m,r]) e of Fig. ||. 
This dressing is achieved without changing the spectral 
asymmetry in any region of Fig. |^ in which the single- 
particle gap remains open, for the U(2) gauge transfor- 
mation is not singular. 

The paramctrization 

b'" = -a,„ 



(4.5) 



sin a cos e ( a^^ - -d^u , , 



-sin a sin 6* ( a^^ - -d^t 



1 - cos a 
cos aa^^ + d^e 



b'll + b'i^Yi'^ where /i = 0, 1, 2 follows from inserting 
Eqs. (|3.4d| ), ( M^ ), and (pTO| ) into Eq. (|^). The trans- 
formation law of Eq. (|4.5 |) under the local U(1)(3)U(1) 



gauge transformation 



cos(2(^5; 
sin(2(/>i 



,>b'^ 



^sm{2<j>,)b'^, 
cos{2<^,)b'l 



(4.6) 



d.. 



At this stage, it is convenient to define the effective 
action (Lagrangian) 



cefi rn' 
m,rj I fj. 



(4.7) 



in the background field B'^ given by Eq. (4.6). This ef- 
fective action is constrained by the gauge symmetries in 
the following way. 

Any transformation of the Grassmann integration vari- 
ables x' and x' with unity for the Jacobi an l eaves the nu- 
merical value of the partition fimction (L2) unchan ged . 
As the Grassmann measure in the partition function (4.2) 
is invariant under the local U(l)(8iU(l) transformation 



Vx!, y := e+*('^-'^=^^3\ (4.8) 



it follows that 

ZLM^ = Z'rn.viV^B'^V - Vhd^V]. (4.9) 
The partition function ([4.2|) thus takes the form 



z:„,,[i?;] = exp 



" '~'m,ri 



(4.10a) 



where 



C. 



off 



^00 



"P*" b"d b" 



b"d b" 



+ C^l^e-'P%°d^b'^ 



(4.10b) 



up to first order in a derivative expansion. This La- 
grangian changes by the usual Abelian Chern-Simons 
boundary terms under the gauge transformation (WM). 



^11 ■ 



and 



The real- valued coefficients C^ii , Cq]^ , ^'33^1 

Cpg' are functions of the parameters m € K. and 77 e M 
with \ri\. A tedious calculation, summarized in Ap- 
pendix yields the values shown in Table |. 



Observe that the coefficients cilt' , ciV , and Ckt.' that 



.(1) 



00 ' ^33 ' 



•'03 



multiply the terms fixed by the local U (1)(8)U (1) gauge 
invariance in the effective Lagrangian ( 4.10b| ) can only 

take a discrete set of values, while the coefficients (7^°' 

and that multiply the terms fixed by the global U(l) 
gauge invariance can vary continuously with m and 77. 

The case m = when TRS is maximally broken is 
special as the symmetry-bre akin g term mSg drops out 
from the Lagrangian in Eq. (4.2V The matrix V in the 
change of Grassmann variables (4.8) is then not restricted 
to the Abelian subgroup U(l)(g) U(l) of U(1)®SU(2) but 
can be arbitrarily chosen in U(2). Consequently, C^]^ = 
C[i' in this limit, which is consistent with the values in 
Table ||. These (equal) coefficients then multiply an SU(2) 
non- Abelian CheritSimons term when m = 0, and hence 
must be quantizcdc3, i.e., 



m— 0,t; 



An 



V p K 



^^e'^P'^b^'d b'° 



(4.11) 
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TABLE L CoefRcients for the effective action in Eq. ( [f.lO| ). The calculation leading to these values is presented in Appendix ^ 



Oil' 



c: 



(1) 



a 



(1) 



C, 



(1) 

03 



\ri\ < m 
m < |?7 



37r|r7| 



4;; sgn77 







where the second hne on the right-hand side is nothing 

but the level 1 SU(2) Chern-Simons term. 

In the case 77 = when TRS holds Eq. ( [4.10bD simpli- 
fies to 



off 

m,r/— 



2tt 



2tt 



e^P^b'^d b'- 



(4.12) 



Notice that the second line is a double Chern-Simons 
term on the fields b'^^ud b'^ which is also called a BF 
Chern-Simons theoryMJlla 

We close this section with the main intermediary step 
of this paper from which the fractionalization of the 
fcrmion charge and s tatis tical phase follows. Insertion 
of Eq. (|4.5|) into Eq. (4.1C) gives the effective action 



= Cj? sin^ a (al ^a^g j (a^, - \dpt 



C« sin^ a e'^P'^ ( a^, - ^-dj ) 9, ( a^, - -dj 



1 



1 - cos a \ / 1 - cos a 
2 '5^^ ) ( cosaag^ + 



1 



c'i^e''P''a^dp ( cos a a 



1 — cos a 



5k 



do 



(4.13) 



with the local U(1)(8)U(1) gauge invariance 



5/^ 



2^, 



(4.14) 



51 



for any compact and boundary- less manifold in (2+1)- 
dimensional space and time. 

Some comments are of order here. First, the coefficient 
controls the axial phase stiffness of the Anderson- 
Higgs contribution to the effective action. Second, each 
of the coefficients Cqq'' , , and C^^^ multiplies a Chern- 
Simons term that is diagonal with respect to the gauge 
fields. The coefficient Cq3'' is different in that regard 
since it couples the gauge field responsible for the 
conservation of the fermion number to the axial gauge 
field flg^ on the one hand, and the axial singlet linear 
combination cig^ = — on the other hand. Such 

an off-diagonal coupling is reminiscent of so-called BF 
Chern-Simons theories.QEJ It is the coefficient that 



controls the charge assignments in the field theory (4.13) 
and. for later convenience, we break its contribution to 



the induced fermionic charge into two pieces. 



BF 

(1) 
BF 

(2) 
BF 



c 



(1) 

BF 



c 



(2) 
BF' 



-"03 

-C^g'adag. 



5i 



(4.15) 



Here, we have introduced the short-hand notation adb ; 



V. 



FRACTIONAL FERMION CHARGE 



Equipped with Eq. (4.13) and Table | we compute in 
this section the leading contributions in the gradient ex- 
pansion to the expectation value of the conserved charge 
current 



if 



Sa,, 



(5.1) 



a, =0 



The induced fermion charge current is 

1 — cos a 



+ 



(5.2) 
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It obeys the continuity equation 
The total induced fermionic charge 



Q := I d'rf{r,t) 



is thus time-independent and given by 

1 — cos a 



Q 



dl ■ cos a Qc 



-d0 



(5.3) 
(5.4) 

(5.5) 



with the help of Stokes' theorem. The induced fermionic 
charge is 



Q = -27r C^3^ f Y + ^ (% - ng) cosa 



(5.6a) 



for the special case when the vector fields and 89 sup- 
port, on a circular boundary at infinity, the net vorticity 



"5 


■ 


2 


2' 






-n 


Al 


1^ ^2 



(5.6b) 



respectively. In the absence of the axial gauge flux rig = 
0, while the condition for the axial vorticity to screen the 



Notice that because C, 



(Kekule) vorticity is = 
vanishes for \ri\ > m, there is no charge bound to the 
vortices in that regime. In contrast, when |?7| < to, the 
charge bound to the topological defect is 



Sm^ f Hg 



-sgn^i, X 



unscreened (rig = 0), 
screened (rig = Ug). 



(5.7) 

These results are consistent with those in Refs. [ij, 
and 



VI. FRACTIONAL STATISTICAL ANGLE 



We start from the effective partition function 



|exp 



" '~'m,Ti 



(6.1) 



with the Lagrangian given by Eq. (4.13) and the coef- 
ficients in Table ||. In a static approximation, i.e., if we 
ignore dyna mics as we did when computing the fractional 
charge (5/7), vortices are independently supported by the 
axial gauge field Og or by the phase 6. 

We will analyze the exchange statistics in two separate 
cases. The first is when the 9 vortices are dynamically 
screened by the half fiuxes in the axial gauge field Og. 
The second case is when the axial gauge field is sup- 
pressed, and the 9 vortex is unscreened; this situation 
does not arise from the effective Lagrangian (4.13) itself. 



but it can occur when one goes beyond the linearized 
Dirac approximation or includes other lattice effects. 



A. Screened vortices 

The exchange statistics of vortices and axial gauge 
fluxes follows from the effective Lagrangian for the so- 
called vortex currents. One way to obtain this effective 
Lagrangian in the screened case is to notice that the lo- 
cal axial gauge invariance together with the first line in 
Eq. (4.13) provides the screening condition, for the axial 
gauge potential must then track the 9 field and, in par- 
ticular, vortices in 9 must be screened by half fluxes in 



One way to impose this screening is to replace 



d^O - 2ag^ 







(6.2) 



in Eq. ( 4.13 ). This can be justified more precisely by 
using the (vortex) dual description of the XY model, as 
presented in Appen dix ^. In effect, the fluctuations away 
from the condition ( |6.2| ), which are penalized by the flnite 

stiffness coefficient j' sin^ a, can be accounted through 
a Maxwell term in the dual description. However, the 
Maxwell term does not enter the exchange statistics. 
Thus, we can simply use t he i nfinite stiffness limit or, 
equivalently, the condition (6.2). 



The La grangian given by Eq. (4.13) in the screening 
limit (U) is 

+ \ Cg^ dp {dJ) (6.3) 



The Lagrangian can be written in terms of the vortex 
current 



that obeys the conservation law 



d V 



0, 



(6.5) 



using the duality representation of the XY model supple- 
mented by a Chern-Simons term in (2-1-1) space and time 
as done in Appendix ^ This leads to the Chern-Simons 
Lagrangian 



-off 



-Cg> {e'-P^d^d^d^ + AndJ:,,) 



1 

4^ 



(6.6) 



from which the statistics carried by screened quasiparti- 
cles with the current j^^.^ follows. This statistics depends 
on the coefficients in Table |. We treat separately the two 
phases of Fig. 
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1. Weak time-reversal symmetry breaking: \r]\ < m 

and the effective La- 



(1) _ 

^00 ~ '-'33 

grangian ( |6.6[) reduces to 



eft 



1 



-sgn/i^a^j^ 



(6.7) 



Thus, because of the absence of the Chern-Simons 
terms, the statistical angle O under exchange of any two 
screened quasiparticles is bosonic, 



TT 



(6.8) 



Notice that it also follows that the induced fermionic U(l) 
current 



1 



-sgn^, Jvrt 



(6.9) 



that couples linearly to a^, is tied to the vortex current. 
In other words, screened quasiparticles with unit vorticity 
are charged objects with charge Q = ±1/2 as found in 
Refs. |lO[ H |l6], and in Sec. 0. 



2. Strong time-reversal symmetry breaking: \ri\ > m 



In this limit, Cw? = Cgg"* = (47r) ^sgnr; and the effec- 



tive Lagrangian (6.6) reduces to 



^sgn77e^'""a^apa^ 



(6.10) 



Using the coefficient of the Chern-Simons Lagrangian for 
the gauge field d^, and its coupling to the vortex cur- 



rent (see appendix ^ for the relation between the 
statistical angle and the coefficient in front of the Chern- 
Simons term), the statistical angle under exchange of 
two screened quasiparticles with unit vorticity is 



e 

TT 



1 



-SgX\T]. 



(6.11) 



Notice that the U(l) current now vanishes, i.e., screened 
quasiparticles carrying fractional statistics are now 
charge neutral. 



B. Unscreened vortices 

We turn to the situation when the axial gauge half 
fluxes are suppressed, while 9 vortices are still present. 
We call these vortices unscreened quasiparticles. This 
situation arises if, in addition to the effective La- 



grangian ( 4.13 ) which followed from integrating out the 
Dirac fermions, there are terms in the effective La- 
grangian due to lattice degrees of freedom that break the 
axial gauge symmetry. For instance, acoustic phonons 
and ripples in graphene can bring about the axial vector 
potential ag ; however, in these cases there is an energy 
penalty of the form Cj^Og that breaks the axial gauge 
invariance due to contributions to the elastic energy. 

The case when the axial gauge potential is absent, i.e., 
the quasiparticles are unscreened, is implemented by the 
replacement 



(6.12) 



in Eq. There follows 



J 



C« sin^ I + cos^ I) sin^l e'^''^dj)d^{dj) 



(6.13) 



(1) 

03 



sm — e 



This Lagrangian can be dualized with the help of the vortex current (6.4) (see appendix 



-C^", = - (Stt^ c[°^ sin' a)-' ^ + J^rt 

•'00 



2'K&^ sin^ I 



(6.14) 



eg' sin^ 



a 



Cii cos^o sin- 
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(The Maxwell term f^J^""^ 
potential c , see appendix 



is associated to the gauge 
s). We shall denote with 



the first line of Eq. (6.14). The statistics carried by un- 



screened quasiparticles with the current j^^.^ follows. This 
statistics depends on the coefficients in Table |[ We treat 
separately the two phases of Fig. 0. 



screened quasiparticles with unit vorticity is 



e 

TT 



_^ sin2a = |^|Q|(l-|g|) (6.16) 
Dm 6m 



by Eq. (5.7). Notice that it also follows that the induced 
fermionic U(l) current, 



Weak time-reversal symmetry breaking: I77I < m 



(6.17) 



In this limit, C, 



(1) 



c: 



(1) 



33 



0, C[l^ = r,/(67rm), C^g 



{2tt) ^sgn /ijj, and the effective Lagrangian ( |6.14 ) reduces 
to 



7] 



247rm 



■ sm a e 



47r d^Jvrt) 



(6.15) 



Using the coefficient of the Chern-Simons Lagrangian for 
the gauge field and its coupling to the vortex current 
j^^t (see appendix ^ for the relation between the statisti- 
cal angle and the coefficient in front of the Chern-Simons 
term), the statistical angle Q under exchange of two un- 



that couples linearly to a^, is tied up to the vortex cur- 
rent. In other words, unscreened quasiparticles with 
unit vorticity are charged objects with charge Q — 
± sin^(a/2) that vari es c ontinuously as a function of the 
ratio jji^/m [see Eq. (^.7|)] as found in Refs. ^ [ij, and 
0. 



2. Strong time-reversal symmetry 



\rj\ > m 



In this limit, C, 



(1) 



— '-'33 — 



(1) 



(4^)-isgn7y, C^;' 



(37?^ - m^)/(127r77^)sgn77, C^^ 
Lagrangian (6.14) reduces to 



= 0, and the effective 



-eff 



sgn 77 



1 - 



377^ 



m 
3i] 



-sm - 



a 



(6.18) 



Using the coefficient of the Chern-Simons Lagrangian for 
the gauge field and its coupling to the vortex current 
(see appendix ^ for the relation between the statisti- 
cal angle and the coefficient in front of the Chern-Simons 
term), the statistical angle under exchange of two un- 
screened quasiparticles with unit vorticity is 



e 

— = sgn 77 

TT 



Sgn 77 



1 - 



3rf 



3jf 



— ^ sm — 
3r;2 2 



sm — 
2 



(6.19) 



\Ql 



Here, we have used the value of the charge \Q\ ~ 
sin^(a/2) for the complementary phase I77I < m. Notice 
that the induced fermionic charge current now van- 
ishes. The unsc reene d quasiparticles carrying the frac- 
tion al sta tistics (S.19) arc thus charge neutral, i.e., \Q\ in 
Eq. ( |6.19 ) should not be confused with the (now vanish- 
ing) electronic charge of unscreened quasiparticles. 

We stress that the quenching of the dynamics in the ax- 
ial gauge field implies the breaking of the axial gauge 
symmetry. It can be thought of as a mean-field approxi- 
mation needed to interpret the numerical simulations of 



the Berry phase acquired by the Slater determinant of 
lattice fermions when one vortex is moved in a quasi- 
static way along a closed curved around another vortex. 
The quench approximation can also be justified if terms 
that break explicitly the axial gauge symmetry such as 



mass term for Og were added to the Lagrangian (4.13) 



After all, from a microscopic point of view, axial gauge 
symmetry is by no means generic. The axial gauge fields 
can be viewed as phonon-induccd fluctuations in the aver- 
age separations between ions that an elastic theory gener- 
ically induces. A mass term for these phonons cannot be 
ruled out by symmetry. 



C. Adding one more fermion to the midgap states 

All calculations for the fractional charge and exchange 
statistics done so far apply at zero chemical potential 
/i = 0, and at some finite staggered chemical potential 
/ig 7^ 0, assuming global vortex neutrality. Global vortex 
neutrality is imposed to bound the energy from above 
in the thermodynamic limit or if periodic boundary con- 
ditions arc imposed. A staggered chemical potential is 
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needed to lift the near degeneracy between the two single- 
particle midgap states that are exponentially localized 
about a vortex and anti-vortex in the bond-density-wave 
(Kekule for graphene) order parameter A, respectively, 
whose separation r is much larger than 1/m. On the 
one hand, when = 0, the two single-particle midgap 
levels are, up to exponentially small corrections in mr, 
pinned to the band center E ~ 0. In the thermodynamic 
limit, their occupancy when /i = is then ambiguous. 
On the other hand, when fi^ ^ 0, the two single-particle 
midgap levels get pushed in opposite directions, one to 
i? > and the other to £^ < (which one goes which way 
depends on the sign of ^g). The single-particle midgap 
level with E' < is then occupied, the other empty, when 
11 = and the results of Sees. ^ and VI B for the frac- 
tional and exchange statistics, respectively, apply. We 
are going to prove that when m > \fi\ > so that the 
two single-particle midgap levels are either both empty or 
both occupied, the exchange statistics is that of semions. 

Suppose one adds one more electron to the Dirac sea 
(here defined to be the Fermi sea at = 0), filling the 
single-particle midgap state at -B > 0. What happens to 
the exchange statistics? 

The easiest way to answer this question is by realiz- 
ing that the Berry phase accumulated by a many-body 
wave function that can be written as a single Slater de- 
terminant (the case in hand) is just the sum of the Berry 
phases for single-particle states. If we fill one more level, 
we only need to add the Berry phase due to that single- 
particle state to that of the filled Dirac sea that we al- 
ready computed. The contribution from the extra level 
can be obtained as follows. (Here we focus on the case 
r] = 0. A generalization to rj ^ can be similarly formu- 
lated.) 

A single-particle midgap wave function is localized near 
a vortex, i.e., its spatial extent is of order 1/m. Details 
on A for distances much larger than 1/m away do not 
matter. Hence, when winding another far-away vortex 
around the first one, the local order parameter A in the 
vicinity of the first vortex just sees its phase change by 
27r. This allows us to focus solely on the problem of 
determining what happens to the single-particle midgap 
wave function as the phase of the order parameter near 
a vortex is rotated by 27r. 

The solution for the single-particle midgap wave func- 
tion when /ig = and in the Dirac approximation was ob- 
tained in Ref. |l^ for the unscreened vortex and in Ref. |l^ 
for the screened vortex. In both cases, the wave function 
picks up a phase of tt when the phase of A changes by 
2tt. If /Xg 7^ 0, the result remains the same, because while 
the midgap level moves with /ig the wave function is in- 
dependent of fig (the wave function has support in only 
one of the sublattices or Ag of the underlying lattice 
model, so the finite value of the staggered chemical po- 
tential does not perturb the single-particle midgap wave 
function). 

In conclusion, occupying one additional single-particle 
fermion level adds a phase of tt to the many-body Berry 
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FIG. 2: (Color online) The induced fermionic charge of a 
quasiparticle, a unit vortex in A(r) with or without attach- 
ment of an axial gauge half flux in a^{r), as a function of 
rj. This charge is computed from the spectral asymmetry of 
spinless fermions hopping on the square lattice with lattice 
spacing o and with a magnetic flux of tt in units of the flux 
quantum (po = hc/e threading each elementary plaquette in 
the static background of a unit vortex in A(r) with or with- 
out attachment of an axial gauge half flux in as{r) and for a 
uniform value of rj. The square lattice is 100 x 100 in units 
of the lattice spacing and the area used for integrating the 
local density of states is a square of size 50 x 50 centered 
around a unit vortex. The following parameters were chosen: 
the hopping t = 1, the magnitude of A(r) on the boundary 
is Ao(oo) = 0.5 while the magnitude of the staggered chemi- 
cal potential is fig = 0.1 (m ~ 0.51). Each (red) filled circle 
is the induced fermionic charge of a unit vortex in A(r) to 
which an axial gauge half flux in 05 (r) is attached. The (red) 
dashed line for rj < m represents the (5=1/2 line. Each 
(blue) filled triangle is the induced fermionic charge of an 
unscreened unit vortex in A(r), i.e., the vector axial flux in 
05 (r) vanishes everywhere. The (blue) solid line for rf < m 
represents Q = 0.402, i.e., the predicted value from the field 
theory with the input parameters. When m ^ rf, the induced 
fractional charge vanishes. A quantum phase transition at 
m = 77, as measured by the jumped in the induced fermionic 
charge, is smeared by finite size effects. 



phase when 77 = 0. This means that the statistical angle 
shifts by A8 = ±7r/2, the statistical angle for a semion, 
when one fermion is added (removed) to (from) the Dirac 
sea. 



VII. NUMERICAL CALCULATION OF THE 
CHARGE AND BERRY PHASE 



We are going to present numerical results on the charge 
and statistics of unscreened vortices supported by the 
bond-density-wave (Kekule for graphene) order param- 
eter A in the presence of the compatible and compet- 
ing order parameters (masses when space and time in- 
dependent) /ig and rf, respectively. The dependence of 
the induced fermionic charge of vortices in A as a func- 
tion of the staggered chemical potential fi^ was studied 
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in Rcfs. |T| and |6| (see also Ref. |^ when =0). The 
foUowing numerical results with the competing mass 77 
are new. 

Our studies have been carried out for the honeycomb 
lattice, which is of direct relevance to graphene, and the 
square lattice with 7r-flux phase. Both lattice models 
yield consistent numerical results. In this paper, only 
the results for tt-Hux phase are presented. The relevant 
technical details for our numerical calculations arc sum- 
marized in Appendix To compare the numerical with 
our analytical results, derived from the Dirac Hamilto- 
nian (ETJ), which is the continuum limit of the linearized 
lattice Hamiltonian, two important issues arise. 

The first one is that all band curvature effects, present 
in any microscopic lattice model, are absent in the con- 
tinuum model. Here we expect that, as long as the char- 
acteristic sizes over which the order parameters vary are 
large compared to the size of the unscreened vortex core, 
static results obtained within the continuum approxima- 
tion should capture some static long-wave length prop- 
erties of the lattice model. This first expectation can 
be concretely addressed by the numerical studies of the 
induced fermionic charge of unscreened vortices below. 

The second issue that arises when one starts from 
a lattice model is the assumed axial gauge invariance 
of Hamiltonian (pTl|). This issue is subtle and sub- 
stantial. The Dirac Hamiltonian (2.1) has a local 
U(l)xU(l) gauge symmetry, while this symmetry is ab- 
sent in graphene, say. Although the vector axial gauge 
field ag is realized in graphene, say through acoustic 
phonons generating ripples, and thus couples in an axial- 
gauge-invariant way to the fermions in the linear approx- 
imation, its kinetic energy is by no means required to be 
gauge invariant. For example, the kinetic energy of is 
expected to contain the axial-gauge-symmctry-breaking 
mass term jag p. It is thus difficult to justify the axial 



gauge invariance of Dirac Hamiltonian ( p.l| ) in a lattice 
model as simple as graphene. 

We do not expect predictions based on Hamilto- 
nian ( |2.l[ ) that rely crucially on the dynamics of the ax- 
ial vector gauge field to capture the corresponding 
low-energy and long wave-length dynamical properties 
of graphene. We will verify this expectation with lat- 
tice computations that require the dynamics of the axial 
gauge field, for example the induced Berry phase as one 
moves a composite particle made of a vortex and an axial 
gauge flux around another composite particle. 

In Sec. VIII, we will present a lattice model that. 



by construction, has the desired local U(l)xU(l) gauge 
symmetry. This model can be used to compute numer- 
ically the statistical phases of unit bond-density-wave 
(Kekule for graphene) vortices screened by axial gauge 
half fluxes and to verify that non-linearities in the many- 
body excitation spectrum do not affect the exchange 
statistics of vortices separated by distances much larger 
than their vortex co re, i .e., this is one model that regu- 
larizes Hamiltonian (pTJ) on the lattice. 



While the system presented in Sec. VIII serves by it 



self as a proof of principle that one can realize the local 
axial gauge invariance on the lattice, the computation of 
the exchange statistics of vortices in this lattice model 
is a computational challenge in lattice gauge theory, as 
opposed to the much simpler exercise in exact diagonal- 
ization for any non-interacting lattice model. 

For this reason, we now limit the numerical studies of 
the statistical phases to the simpler case when Og ^ 0, 
i.e., the case of unscreened vortices. In effect, we are 
ignoring all many-body effects imposed by the local ax- 
ial gauge invariance and thus treating the problem at the 
mean-field level. By comparing the charge obtained from 
the Aharonov-Bohm effect with that obtained directly 
from the local density of states, we will show that this 
approximation is qualitatively (but not quantitatively) 
justified for dynamical properties of bond-density-wave 
(Kekule for graphene) vortices, whereas it fails dramat- 
ically for dynamical properties of the axial gauge half 
fluxes. 



A. Static calculation of the charge 

We begin with the study of static properties, when the 
vortices or axial gauge half flux tubes are not moved, 
so that the dynamics of the axial gauge potential is not 
relevant. One physical quantity that can be studied in 
the static limit is the induced fermionic fractional charge. 
It is obtained by summing up the local fermionic density 
of states in a region of space that encloses the core of the 
vortex. 

In our numerical studies, a vortex is placed at the cen- 
ter of the square lattice system of size 100 x 100 in units 
of the lattice spacing a while a flux of n in units of the 
flux quantum (/)g = hc/e threads each elementary plaque- 
tte. An area of integration, 50 x 50, centered around the 
vortex is used for summing the local fermionic density of 
states. We fixed the strength of the bond-density-wave 
(Kekule for graphene) order parameter A = 0.5 and stag- 
gered chemical potential fi^ =0.1. 

Figure || shows the value of the induced fermionic 
charge as a function of rj with and without the axial 
gauge half fiux. A clear normalization of the fractional 
charge to 1/2 follows from adding an axial gauge half 
flux. Notice that there is a (smoothed) step as the mass 
77 becomes comparable to m. This is the finite-size sig- 
nature of a quantum phase transition at \r]\ = m. The 
results in Sec. are displayed in Fig. |[ They correspond 
to sharp step functions at the transition point I77I = to. 
The numerical results displayed in Fig. |^ are consistent 
with the analytical results ( |).7[ ) , keeping in mind that the 
lattices studied are finite and thus quantum transitions 
are smeared. For that matter, notice that the agreement 
between the field-theory prediction and numerics is best 
away from the critical point I77I = m. 
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FIG. 3: (Color online) Schematics of the static magnetic 
flux tubes inserted to probe the induced fermionic charge of a 
quasiparticle, a mass vortex with or without the attachment 
of an axial gauge half flux, using the Aharonov-Bohm effect 
in the second set-up described in the text, (a) We insert 
one static magnetic flux tube (colored in red) with the flux 
(j) = l(j>o (the flux quantum is (/jq — hc/e) while a quasiparticle 
encircles dynamically this magnetic flux with the trajectory 
indicated by the directed loop (colored in blue) . (b) We insert 
two static magnetic flux tubes (colored in red) with the fluxes 
cj) = ±l<j>o while a quasiparticle encircles dynamically one and 
only one magnetic flux tube with the trajectory indicated by 
the directed loop (colored in blue). 



B. Dynamic calculation of the charge 

A dynamical alternative to computing the induced 
fermionic charge through the integrated local density of 
states is the following. If we take a unit vortex in A (with 
or without an accompanying axial gauge half flux) around 
a circle of radius r that encircles a magnetic flux, then 
an Aharonov-Bohm phase accumulates. The value of the 
charge induced near the vortex follows after matching 
the Berry phase computed numerically to the analytical 
value of the Aharonov phase. 

We carry out this approach in two different set-ups. In 
the first, we apply a uniform magnetic field to the system, 
i.e., we fix a given electromagnetic flux 



^00 > 



Z € R, 00 the quantum of flux, 



(7.1) 



per elementary unit cell on the lattice. The Aharonov- 
Bohm phase 7j^g that is picked up depends on the radius 
of the path since the encircled magnetic flux scales with 
the area. The Aharonov-Bohm phase in this case is thus 
given by 



-y^g = 2tt X Q X {nr'^) x I. 



(7.2) 



Here, Q is the charge bound to the unit vortex in A. 

A second set-up is shown in Fig. ^a. We insert an 
electromagnetic flux tube with flux 



through the elementary unit cell on the lattice at the cen- 
ter of the system. All other elementary unit cells are free 
of any magnetic flux. We then move the unit vortex in the 
bond-density-wave (Kekule for graphene) order parame- 
ter around a path enclosing this flux. Notice that the 
Aharonov-Bohm phase is independent of the path 
as long as it strictly contains the magnetic flux tube, i.e., 
the elementary unit cell at the center of the lattice. It is 
expected to have the value 



7ab = 271 X Q X I. 



(7.4) 



We also study the case displayed in Fig. ^jb. The rea- 
son for it is that we want to ensure that compensating 
fermionic charges on the edges of the sample do not con- 
tribute a phase as well. In the set up of Fig. ||b, whatever 
happens with the fermionic edge charges does not lead 
to an Aharonov phase because their path would encircle 
(even if they move) the vanishing total flux 



(7.5) 



The results we obtain for the Berry phase when we 
wind the unscreened vortices around a closed path are 
shown in Fig. |^ for the case of the first set-up (uniform 
applied magnetic field). We fix the parameters A = 0.5, 
r = 14.5 (in a 56 x 56 lattice) and </> = 0.001(/)q per pla- 
quette, and plot the charge Q versus the parameter /x^/A. 
The blue dots and red dots are the numerical results for a 
vortex without the axial gauge half flux and with the ax- 
ial gauge half flux, respectively, while the corresponding 
theoretical predictions from Ref. and 16 are plotted 



in blue and red solid line. Notice that the analytical and 
numerical results agree quite well for the case of vortices 
unscreened by axial gauge half fluxes. 

As anticipated; the analytical and numerical results 
are not consistent for the case of screened vortices. The 
reason is precisely what we highlighted in the beginning 
of this section, i.e., that the lattice model studied nu- 
merically in this section does not contain the U(l)xU(l) 
symmetry, i.e., the axial gauge fleld dynamics present 
in the Dirac Hamiltonian (2.1). The same issue applies 
to the problem of computing the exchange statistics of 
pairs of screened vortices. We cannot study the statis- 
tical angle of screened vortices within the approach of 
this section. In Sec. VIII, we will present a microscopic 



model that does have the U(l)xU(l) gauge symmetry. 
However, this model cannot be studied by simply com- 
puting Slater determinants (see Appendix ^ as has been 
done so far in this section. 



^00 > 



I G 



the quantum of flux. 



(7.3) 



Before closing Sec. VII B, let us mention that we have 
checked the results summarized by Fig. ^ that we ob- 
tained by applying a uniform magnetic field against those 
obtained with a single flux tube as in Fig. ^ or with two 
flux tubes as in Fig. ||b. 
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FIG. 4: (Color online) The induced fermionic charge of a 
quasiparticle, a unit vortex in A(r) with or without attach- 
ment of an axial gauge half flux in a^{r), as a function of the 
ratio /ig/Ao(oo). This charge is obtained by matching the 
numerical Berry phase picked up when a quasiparticle hops 
along the closed boundary of an area that encloses a magnetic 
flux in the uniform background of (j.^ to the correspo nding 
Aharonov-Bohm phase along the lines outlined in Sec. VII B 
and Appendix |^. Hopping takes place on the square lattice 
with the lattice spacing a and a magnetic flux of tt in units 
of the flux quantum (j>o = hc/e per elementary plaquette, 
there being 56 x 56 elementary plaquettes. The closed path 
used to compute the Berry phase is approximately circular 
with the radius r = 14.5 in units of the lattice spacing. The 
following parameters were chosen: the hopping t — 1, the 
magnitude of A{r,t) on the boundary is Ao(oo) = 0.5 while 
the flux is (p = 0.00100 through each elementary plaquette. 
The (red) flUed circles are the induced charges of a dynamical 
unit vortex in A{r,t) to which is also attached a dynamical 
axial gauge half flux in a^{r,t) as a function of ^s/Ao(oo). 
The (red) dashed line is the analytical charge Q — 1/2. The 
(blue) flUed triangles are the induced charges of a dynami- 
cal unit vortex in A(r,t) without the axial gauge half flux 
in a5(r) as a function of the ratio ^s/Ao(oo). T he ( blue) 
solid line is the induced charge computed from Eq. (5.7) as a 
function of the ratio fi^^/ Aa{cx>). 



C. Fractional statistics for unscreened vortices 

We now present the numerical value of the statisti- 
cal angle Q in units of tt acquired under the exchange 
of two unit unscreened vortices in the bond-density-wave 
(Kckulc for graphcnc) order parameter A, which we shall 
call quasiparticles from now on. Wc have computed nu- 
merically the Berry phase 7 in units of tt accumulated 
when a first dynamical quasiparticle moves along a tra- 
jectory that winds once around a second static quasipar- 
ticle as outlined in Appendix The statistical angle O 
acquired under the exchange between these two quasi- 
particles is then 



pcndix P and unlike in Sec. VII l| , we only treat un- 
screened vortices. We have verified that 7, when com- 
puted along the lines of Appendix ^ docs not change 
when axial gauge half flux tubes are attached to the vor- 
tices. 

To compare the microscopic exchange statis tics w ith 
the one computed within field theory in Sec. VI E| , we 
restrict the numerical computation to the half-filled case. 
However, we will also test the prediction of Sec. VIC 
by working with one spinless fcrmion more than (or less 
than) at half-filling. 

We will always take m w 0.5 in Eq. ( |2.2| ). Moreover, 
to limit finite size effects, we assume that 77 ^ to, i.e., we 
work well below the transition point jryj = to when the 
breaking of TRS is weak. 

The 77 dependence of the Berry phase 7 with to fixed is 
shown in Fig. |^ for different values of the uniform stag- 
gered chemical potentials /ig. The magnitude of the Berry 
phase is seen to be independent of whether the pair of 
quasiparticles have the same (filled circles) or opposite 
(star symbols) vorticities, but it does depend on /i^, i.e ., 
on the induced fractional charge Q given in Eq. (^.7|) . 
The r] dependence of 7 is linear, as predicted in Sec. |VI B| , 
but with slope s dev iating from the theoretical predic- 
tions, i.e., Eq. ( 3.16 ), shown as the solid or dashed lines. 
The agreement betwe en th e Berry phase of the micro- 
scopic model and Eq. (6.16) is thus qualitatively but not 
quantitatively good. 

The microscopic Berry phase 7 as a function of the 
ratio A/to, which also parametrizes (3(A, ^g), when ry = 
0.025 is held fixed is shown in Fig. || as filled circles when 
the quasiparticles carry the same vorticities or as stars 
when the quasiparticles carry the opposite vorticities. As 
expected, exchanging a pair of quasiparticles with equal 
unit vorticities differs solely by a sign relative to exchang- 
ing a pair of quasiparticles with opposite unit vorticities. 
The lines (solid when the quasiparticles have the sam e 
unit vorticity, dashed otherwise) are given by Eq. (3.16). 
Evidently, the dependence on Q ^ 1/2 of the microscopic 
exchange statistics is not captured by the field theory. 



(7.6) 



Here, as we do not impose dynamically the axial gauge 
symmetry at the microscopic level as presented in Ap- 



As discussed in Sec. VI C , when adding (removing) one 
fcrmion to (from) half-filling, the Berry phase accumu- 
lated by a complete winding of quasiparticles of opposite 
unit vorticities changes by tt for the case 77 = 0. This 
extra phase is the response of the single-particle midgap 
states to varying the phase of A by 2tt. Numerically, this 
assertion is confirmed directly by computing the accumu- 
lated Berry phase and obtaining 7 = ±7r when filling or 
emptying one midgap state. 

In summary, comparison of the microscopic Berry 
phase accumulated by winding an unscreened quasipar- 
ticle around a static one with the field -theo ry computa- 
tion of the exchange statistic in Sec. VI B| shows that: 
1) The microscopic Berry phase 7 (and consequently the 
microscopic exchange statistical angle 8 — 7/2) varies 
continuously as a function of rj and in a linear fashion 
for small 77, in good agreement with the field-theory re- 
sults. 2) The slope 7/77 shows a monotonic dependence on 
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FIG. 5: (Color online) Berry phase in units of tt as a function 
of ri/m <^ 1 for fixed m acquired during the exchange of two 
unscreened quasiparticles, i.e., unit vortices in A(r, t) without 
the attachment of axial gauge half fluxes in a^(r,t ). Num er- 
ical computations along the lines outlined in Sec. VII B and 



Appendix ^ were performed for spinless fermions hopping on 
the square lattice with lattice spacing a and with a magnetic 
flux of TT in units of the flux quantum (j>o = hc/e thread- 
ing each elementary plaquette in the dynamic background of 
a unit vortex in A{r,t) without the attachment of an axial 
gauge half flux in a^{r,t) and for a uniform value of rj. The 
square lattice is 72 x 72 and the exchange path is approxi- 
mately circular with the radius r = 18.5 in units of the lattice 
spacing. The following pa rameters were chosen: the hopping 
t = 1, m = Aq(c>o) -f ~ 0.51 but with two different value 
of =0.1 and 0.025. Filled circles and solid lines represent 
the case when the two quasiparticles carry the same unit vor- 
ticity. Stars and dashed lines represent the case when the two 
quasiparticles carry the opposite unit vorticity. Symbols are 
obta ined numerically while the lines are the predictions from 
Sec. IvIBI. 
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FIG. 6: (Color online) Berry phase in units of tt as a func- 
tion of Ao(oo)/m for fixed m and r\ <^m acquired during the 
exchange of two unscreened quasiparticles, i.e., unit vortices 
in A(r,t) without the attachment of axial gauge half fiuxes 
in ac.(r . t'). N umerical computations along the lines outlined 
in Sec. VII B and Appendix ^ were performed for spinless 
fermions hopping on the square lattice with lattice spacing 
a and with a magnetic flux of tt in units of the flux quan- 
tum ^0 = hcje threading each elementary plaquette in the 
dynamic background of a unit vortex in A{r,t) without the 
attachment of an axial gauge half flux in a^{r,t) and for a 
uniform value of rj. The square lattice is 72 x 72 and the ex- 
change path is approximately circular with the radius r = 18.5 
in units of the lattice spacing. The f ollowing para meters were 
chosen: the hopping t = 1, m = \J Aq(oo) + /i^ = 0.51 and 
77 = 0.025. Filled circles and solid lines represent the case 
when the two quasiparticles carry the same unit vorticity. 
Stars and dashed lines represent the case when the two quasi- 
particles carry the opposite unit vorticity. Symbols are ob- 
tain ed nu mericallv while the lines are the predictions from 
Sec. IVIBI. 



the ratio A/m, which is not in good quantitative agree- 
ment with the field-theory results. 3) The magnitude I7I 
is independent of the relative sign of the quasiparticles 
vorticities. This is expected for a vortex and its anti- 
vortex can annihilate. Consequently, winding a third 
vortex around a vortex anti- vortex pair must accumulate 
a vanishing Berry phase. 4) Microscopic semion statis- 
tics O = ±7r/2 is obtained when adding (removing) one 
fermion to (from) the half-filled system in agreement with 
the prediction from the continuum theory. 



VIII. MICROSCOPIC MODEL 



ther on the honeycomb or 7r-flux lattices used in Sec. VI] 



We have seen in Sec. VII that the fractional charge 



induced by an axial gauge half flux in cannot be mea- 
sured dynamically from the Aharonov-Bohm phase in- 
ferred from the numerical computation of a Berry phase. 
This is so because the local axial gauge symmetry in the 
continuum Hamiltonian (2.1) is not present in the lattice 
model used in Sec. VII. Thus, there is a dynamical con- 



tribution that is missing and that cannot be captured by 
the simple models of one species of fermions hopping ei- 



For the same reason, we could not obtain numerically 
the exchange statistics in the case when the vortices are 
screened by the axial gauge potential, since the exchange 
of the topological defects necessarily acquires a dynami- 
cal contribution from Og. 

We now construct a lattice model sharing the same 
local U(l)xU(l) symmetr y an d the same particle content 
as the dynamical theory (3T). The predictions fo r the 
exchange statistics of screened vortices done in Sec. VIA 



should be captured by this lattice model. Unfortunately, 
we cannot verify this claim, for the largest system sizes 
that we could treat numerically are of the order of the 
vortex core. 

Consider a square lattice A whose sites we denote 
with the Latin letters i,j,k, and I. We denote with 
1 = X, 2 = y the two orthonormal vectors spanning the 
square lattice A (and we will index these two vectors as 
fi — 1,2, for fi — 1,2). Links (or bonds) on the square 
lattice between nearest-neighbor sites i and j are labeled 
by (ij) (or simply by ij when used as an index to a field 
defined on the links). We denote by 0-jf.i the square 
plaquette with the corners i, j, k, and I. 

We define four sets of operators. There are the bosonic 
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operators A^^ and A^^^j living on the links of the square 
lattice A. There are the bosonic operators 0^ and the 
fermionic operators ipj_ living on the sites. The spinor- 
valued operator ip^ has here four components on which 
the 4x4 matrices defined in Eqs. (2. Id) and (2.1e) act. 

These four sets of operators, together with their canon- 
ical conjugate operators, satisfy the following relations: 



(8.1a) 



(8.1b) 



--*<S,„-, (8.1c) 



and, finally, 



(8.1d) 

with the equal-time global constraint (half-filling con- 
straint) 



.Ic) 



(Since we are working with four flavors of fermions, half- 
filling means average 2 particles per each site.) 

We define the lattice model by the quantum Hamilto- 
nian 



H ■.= Hg + Hg^+Hj + H,+Ht,+H,^. (8.2a) 



Here, 



2 ' o 



(8.2b) 

describes a U(l) lattice gauge theory with gauge coupling 



9i 



(8.2c) 

describes another U(l) lattice gauge theory with gauge 
couphng .gf , 



1 



-»(<^.-<^,)+2jA,.^ 4,H.c 



.2d) 



describes a quantum rotor (XY) model with coupling J^, 
and 



ieA 11=1,2 
+ H.C. 



.2e) 



describes the nearest-neighbor hopping with the real- 
valued amplitude t of 4 independent fermions per site. 
So far, there are 4 non-equivalent Dirac points at half- 
filhng which are located at fc = (0,0), (0,7r), (tt, 0), and 
(tt, tt). This is why we have added the term 



i?e*^'('+A)+*^5^5,(,+,) ,/,(^^.^ + H.c. 



.2f) 



that opens a gap of order t' at the points k = 
(0, tt), (tt, 0), (tt, tt), thus leaving k = (0,0) as the sole 
Dirac point. This scheme is precisely Wilson's procedure 
used to piiercome the doubling problem in lattice gauge 
theories.Lj An important comment is in order, however. 
One reason why this prescription is not fully satisfying 
in lattice gauge theories is that a ny m ismatch between 
the first and second terms of Eq. ( ) leads to a gap at 
k = (0,0) as well, i.e., fine-tuning is needed to achieve 
the correct particle content. Here, this is fine because we 
are interested in systems where there is such a gap. No- 
tice in that regard that the gap at fc = (0, 0) that arises 
from a small mismatch between these two terms (a small 
fraction of t') is much smaller than the one at the edges 
of the Brillouin zone (order t'). Indeed, such a term due 
to a mismatch is actually part of the final term that we 
consider in the Hamiltonian, namely 

^™ ■■= E ^1 (^s^ + Ak/3e*''^^- + 1^10^2) k ■ (8.2g) 

This contribution docs indeed open a gap at the remain- 
ing Dirac point at fc = (0, 0). 

For any smooth and static boson background, the con- 
tinuum limit of Hamiltonian (8.2) upon linearization of 



the fermion spectrum at the two non-cquivalcnt Dirac 
points at half-filling is given by Eq. (2.1), as is also the 
case with the fermion spectrum of graphene restricted to 
spinless fermions hopping with sufficiently smooth mod- 
ulations of the hopping amplitudes. 

Co ntra ry to graphene for spinless fermions, Hamilto- 
nian (8.2) is invariant under the local U(l)xU(l) gauge 



transformation 



^1 



5y 



A^^ > A, 



A, 



tj {Xi Xj) 1 



(8.3a) 
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generated by 



with 



..3b) 



ieA L V ^=1,2 
( iWs^z + 2n, + ^ L 



i(i+fl) I Xj 



5i(i+/i) I 



(8.3c) 

where Xi s-nd are arbitrary real-vahied numbers. 

The physical subspace is the set of gauge invariant 
states, i.e., states that are tensor products of states in 
the Fock space generated by the algebra Eqs. (O) , 



1^, 



such that Gauss law holds globally, 

G-\x,o\^) = m 



Aa) 



(8.4b) 



for all real-valued function x ^, or, equivalently, lo- 
cally 







+V'h5V', +2nj |vi/), 



(8.4c) 



for any z g A. 

We denote by a gauge invariant state (B.4) with 

two fractional charges localized around sites i and j, re- 
spectively. The statistical phase 6 induced by the physi- 
cal process by which two fractional charges are exchanged 



is given by the difference between two Berry phases,E 



e:=iarg J] 

1 I 
- 2arg n 



(8.5) 



i, 



For both Berry phases, one fractional charge hops along 
the closed path V = {i^}, while the other fractional 
charge is static. For the former Berry phase, j is lo- 
cated inside the area bounded by a choice that we 
denote by j C V . For the latter Berry phase, j is located 
outside the area bounded by 7^, a choice that we denote 
by 2 C V. 

The dimensionality of the gauge-invariant Hilbert 
space scales w ith th e dimensionality of the fermionic 
Hilbert space (8. Id), which itself scales exponentially 



fast wit h the number of sites. Given the half-filling con- 
straint (B.lc), this limits the numerical evaluation of the 
right-hand side of ( |8.5| ) to lattices with linear dimen- 
sions of the order of the core size Xjm of the defects, 
i.e., on dist ances much too short for the right-hand side 
of Eq. (8.5) to be interpreted as the statistical angle of 
point-like quasiparticles. 

If we are will ing to give up the local U(l)xU(l) gauge 
the strongly correlated nature of 



3), i.e. 



mvariance 

the problem, we can compute the contribution to the 
statistical phase arising from the fermion hopping. In- 
deed, the problem then reduces to a single-particle one for 
which the dimensionality of the relevant Hilbert spaces 
only scales linearly with the number of sites. We stress 
that this contribution alone violates the local U(l)xU(l) 
gauge invariance. 



J 



IX. MORE SPECIES OF FERMIONS - 
CLASSIFICATION OF ALL MASSES IN 
GRAPHENE AND tt-FLUX PHASE 

So far we have ignored the spin- 1/2 quantum number 
of electrons. If so, in the linear approximation (2.1) of 

graphene restricted to spinless fermions say, ^^gcaiar ex- 
hausts all possible symmetry-breaking instabilities with 
a local order parameters compatible with charge conser- 
vation. The local order parameter for a charge-density 
wave that breaks the sublattice symmetry but preserves 
the time-reversal symmetry is the real-valued order pa- 
rameter fis{r) (introduced by Semenoff for graphene in 
Ref. ^. The local order parameter for a bond-density 
wave instability that preserves the sublattice and time- 



reversal symmetries is the complex-valued order param- 
eter A(r) (the U(l) Kekule order parameter introduced 
by Hou et al. for graphene in Ref. |l^. The local order 
parameter for a bond-density wave instability that breaks 
the sublattice and time-reversal symmetries is the real- 
valued order parameter r]{r) (introduced by Haldane for 
graphene in Ref. ||). 

If we reinstate spin- 1/2 in the most naive way and con- 
sider two independent copies of the model in Eq. ( [2.l| ), 
then the results we found for spinless electrons are mod- 
ified in a trivial way. Defects bind equal values for 
the fractional charge for both species, up and down 
spin, thereby doubling the total induced fermionic charge 
(which is to be associated with a spin-singlet state). The 
same happens to the exchange statistical angle. It is sim- 
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TABLE II : Th e 36 mass matric es w ith particle-hole symmetry (PHS), see Eq. (9.6), for the massless Dirac Hamiltonian /Cq 
from Eq. (9.4b) are of the form (9.3) and anticommute with /Cq. Each mass matrix can be assigned an order parameter for 
the underlying microscopic model, here graphene or the square lattice with 7r-flux phase. The latin subindex of the order 
parameter's name corresponds to the preferred quantization axis in SU(2) spin space. The pair of numeral subindices 02 and 32 
are used to distinguish the two unit vectors spanning two-dimensional space. Each mass matrix preserves or breaks time-rev ersa l 
symmetry (TRS), see Eq. (3.7), spin-rotation symmetry (SRS), see Eq. (3. J), and sublattice symmetry (SLS), see Eq. (19. 91) . 
To any of the 36 mass matrices corresponds a "partner" mass matrix obtained through the involutive transformation ( |9.12| ) 
denoted C. 



Mass matrix 


Order parameter 


TRS 


SRS 


SLS 


Partner by C 


Order parameter by C 


C invariant 


-^^3010 


ReVBS 


True 


True 


True 


V3010 


ReVBS 


True 


-^^0020 


ImVBS 


True 


True 


True 


V0020 


ImVBS 


True 


-^^3033 


CDW 


True 


True 


False 


V3333 


Neel, 


False 


-^^3003 


QHE 


False 


True 


False 


V3003 


QHE 


True 


-^^3110 


ReVBS, 


False 


False 


True 


V2132 


ImTSC32. 


False 


Y 

^0210 


ReVBS 


False 


False 


True 


Y 

^1132 


ReTSC^o 


False 


^3310 


ReVBS, 


False 


False 


True 


V3310 


ReVBS, 


True 


Y 

^0120 


ImVBS^. 


False 


False 


True 


Y 

^1102 


ReTSCno 


False 


-^^3220 


ImVBS^ 


False 


False 


True 


V2102 


ImTSCo2, 


False 


-^^0320 


ImVBS, 


False 


False 


True 


V0320 


ImVBS, 


True 


-^^3103 


QSHEa; 


True 


False 


False 


^"2121 


ImTSC, 


False 


-^^0203 


QSHEy 


True 


False 


False 


V1121 


ReTSC, 


False 


-^^3303 


QSHE, 


True 


False 


False 


V3303 


QSHE, 


True 


-^^3133 


Neel, 


False 


False 


False 


V2211 


ReSSC 


False 


-^^0233 


Neelj, 


False 


False 


False 


V1211 


ImSSC 


False 


-^^3333 


Neel, 


False 


False 


False 


V3033 


CDW 


False 


-^^2211 


ReSSC 


True 


True 


False 


V3133 


Neel, 


False 


-^^1211 


ImSSC 


False 


True 


False 


^"0233 


Neely 


False 


-^^1002 


ReTSCo2i/ 


True 


False 


True 


^1002 


ReTSCo2y 


True 


-^^2002 


ImTSCoaj, 


False 


False 


True 


V2302 


ImTSCo2. 


False 


-^1102 


ReTSCoa. 


False 


False 


True 


V0120 


ImVBS, 


False 


-''^2102 


ImTSCoa, 


True 


False 


True 


V3220 


ImVBSy 


False 


-^^1302 


ReTSCo22; 


False 


False 


True 


V1302 


ReTSCo2a; 


True 


-^^2302 


ImTSCo2. 


True 


False 


True 


V2002 


ImTSCo2« 


False 


-^^1032 


ReTSC32y 


False 


False 


True 


V1332 


ReTSC32a; 


False 


-^^2032 


ImTSCga^ 


True 


False 


True 


V2032 


ImTSC32j, 


True 


-^^1132 


ReTSCga. 


True 


False 


True 


^0210 


ReVBSy 


False 


-^^2132 


ImTSCga, 


False 


False 


True 


V3110 


ReVBS, 


False 


-^^1332 


R,eTSC32:E 


True 


False 


True 


V1032 


ReTSC32y 


False 


^2332 


ImTSC32. 


False 


False 


True 


V2332 


ImTSC32. 


True 


-^^1021 


ReTSCy 


True 


False 


False 


V1321 


ReTSC, 


False 


-^2021 


ImTSCy 


False 


False 


False 


V2021 


ImTSCy 


True 


-^^1121 


ReTSC, 


False 


False 


False 


V0203 


QSHEj, 


False 


-^^2121 


ImTSC, 


True 


False 


False 


V3103 


QSHE, 


False 


-^^1321 


RcTSC, 


False 


False 


False 


V1021 


ReTSCj, 


False 


-^^2321 


ImTSC, 


True 


False 


False 


V2321 


ImTSC, 


True 



ply doubled with respect to the results in Sec. VI 



However, if spin is not a good quantum number, a 
larger number of instabilities can occur and more masses 
or order parameters (other than Re A, ImA, fi^, and ?/) 



need to be taken into account. Thus, one must consider 
more generic Dirac Hamiltonians and study all their al- 
lowed masses. Topological defects in these order param- 
eters could bind states, whose (fractional) charge and 
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TABLE III: Enumeration of the 56 distinct 5-tupIets of maximally pairwise anticommuting PHS X^^_^^^^^^^^^. 



are broken into 28 pairs related by the operation of C conjugation (9.15 



The 56 5-tuplets 



5-tuplet 

{ReVBS, ImVBS, ReSSC, ImSSC, CDW} 

{imVBS, CDW, ReVBS,, Re VBSy, ReVBS, 
(ReVBS, CDW, ImVBS,, ImVBS^, ImVBS, 
{ReSSC, ImSSC, QSHE,, QSHE^^, QSHEj 



Partner 5-tuplet by C conjugation 

{ReVBS, ImVBS, Neel,, Neel^, Neel^ 

{ImVBS, Neel, , ImTSCga, , ReTSCaa, , ReVBS, } 
{ReVBS, Neel, , ReTSCoa, , ImTSCoa, , ImVBS, } 
{Neel, , Neely , ImTSC, , ReTSC, , QSHE, } 



32z 
32z i 



[ReVBS, ReSSC, ReTSCoa,, ImTSCozy, ReTSC, 
ReVBS , ImSSC, ImTSCoa, , ReTSCoa^ , ImTSC, 
ImVBS, ImSSC, ReTSCga, , ImTSCga^ , ReTSC, 
ImVBS, ReSSC, ImTSCga, , ReTSCaa^ , ImTSC, 
CDW, ImSSC, ImTSC,, ReTSCy, ImTSC, 

[cDW, ReSSC, ReTSC, , ImTSC^,, ReTSC, 



{ImVBS, , QSHEj, , ImVBS, , ReTSCga^ , ImTSCsa^ } 
{ImVBS, , QSHEj^ , ReVBS, , Neel, , QSHE, } 
ilmVBS, , ReTSCsay , ImTSCaa, , ImTSCoa, , ImTSC, 
jlmVBS, , ReTSCga, , ReTSCoa, , ReTSC, , ImTSCgzy 
{ImVBS, , ReTSCsa, , ImTSCga, , ImVBS^ , QSHE, } 
{ImVBS,, ReTSC,, ImTSC,, CDW, ReVBS,} 

{QSHEj^ , ImVBS, , QSHE^ , ReVBS, , Neel, } 
{QSHEj,, ReTSCo2a, ReTSCy, ImSSC, ImTSCgzy } 
{QSHEj,, ReTSCo2a, ImTSCo2j;, ReVBS,, ReVBS, } 
{QSHEj,, ReTSCgay, ImTSCoay, ImTSCj,, ReSSC} 

(ReVBSj, , Neely , QSHE, , ImVBS,, , QSHE _ } 
{ReVBSy, ReTSCy, ImTSCy, CDW, ImVBS J 
{ReVBSy , ReTSCgay , ImTSC^ , ImTSCo2, , ImTSCoa, } 
{ReVBSy , ReTSCoa, , ImTSCoa, , QSHE, , ReVBS, } 

Neely , ReTSCgay , ImTSCoa^ , ImTSC, , ImTSC, } 
ImVBS, , ReTSCsa,; , ImTSCoa, , ImTSC, , ImTSCga, 
ReTSCoay , ReTSC^ , ImTSCga, , ImTSCga, , ImVBSy 
[ ReTSCy , ReTSCoa, , ImTSC, , ImTSCga, , Neel, } 



{ReVBS, Neel, , ReTSCoa, , ImTSCoa, , ImVBS, } 
{ReVBS, Neel^ , ImTSCoa^ , ReTSCoa^ , ImVBS^ } 
{ImVBS, Neel^, ReTSCgay, ImTSCgay, ReVBS J 
{ImVBS, Neel, , ImTSCga, , ReTSCga, , ReVBS, } 
{Neel, , Neely , ImTSC, , ReTSC, , QSHE, } 
{Neel, , Neel, , ReTSC^ , ImTSC^ , QSHE^ } 

{ReTSCoa, , ReTSC, , ImVBS, , ReTSCga, , ImTSCga^ } 
{ReTSCoa, , ReTSC, , ImTSCga, , ReSSC, QSHE, } 
{ReTSCoa, , ReTSCga, , ReVBS, , ImTSCoay , ImTSC, } 
{ReTSCoa, , ReVBS^ , ReTSCoa, , ReTSC^ , ImTSCga^ } 
{ReTSCoa, , ReVBS^ , ReVBS, , ImTSCoa, , QSHE, } 
{ReTSCoa, , ReTSC^ , ImTSC, , Neel, , ImTSCga, } 

{ReTSC, , ImVBS, , ImTSC, , ReVBS, , CDW} 
{ReTSC, , ReTSCoay , ReTSC, , Neel^ , ImTSCga^ } 
{ReTSC, , ReTSCoay , ImTSCoa, , ImTSCga, , ReVBS, } 
{ReTSC, , ReTSCga, , ImTSCoa, , ImTSC^ , Neel, } 



{ReTSCga, , ImSSC, ImTSC, , ImTSCoa, , QSHE, } 
{ ReTSCga, , ReTSC, , ImTSC^ , Neel, , ImTSCoa, } 
I ReTSCga, , ReTSCga, , ImTSC^ , ImVBSy , ImTSCoa 
{ReTSCga, , ReTSCoa, , ImTSCoa^ , ImTSC, , ReVBS 



{ImSSC, ReTSCga,, ImTSCoa,, QSHE^, ImTSC,} 
llmVBS,, ReTSCga,, ImVBSy, QSHE^, ImTSCga,} 
{ ReTSCoaa , ReTSC, , ReVBS, , ImTSCga, , ImTSCoa, } 
{ReTSC,, ReTSCoa,, QSHE,, ImTSCga,, ReSSC} 



statistics would depend on the effective action (as func- 
tion of all the mass order parameters and the a„ and 



microscopic symmetries. 

This classification applies as well to the microscopic 



fields) that is obtained upon integrating all the species of 
fermions. This effective action would be the extension of 
the one derived in Sec. [V for the case of the four order 
parameters (ReA, ImA, /ig, and 77). 

We do not fully carry this program in this paper. 
Nonetheless, we classify all these masses according to the Wilson masses can also classified as we do below. 

I 



model of Sec. VIII. There, we c hose a specific way to 
add Wilson masses [see Eq. (8.2f)] to selectively get rid 
of all but 2 Dirac points in ord er t o recover in the long- 
wavelength limit Hamiltonian (2.1). The set of all (64) 



Classification of masses in graphene and tt-Hux 
phases 



To describe all symmetry-breaking instabilities with a 
local order parameter in graphene or the square lattice 
with TT-flux phase, we consider the Bogoliubov-de Gennes 



(BdG) Hamiltonian 



where *]/ is the 16-component Nambu spinor 



■= (■'/'t, V'j., i'], ■'/'x) 



(9.1a) 



(9.1b) 



21 



and ■i/'s=T i. ^ 4-componeiit fermion annihilation oper- 
ator that accounts for the 2 valley and the 2 sublattice 
degrees of freedom. The kernel of the BdG Hamiltonian 
has the block structure 



jr — I ^PP ^ph 

^ ' '•ph '■•pp 



(9.1c) 



where the 8x8 blocks Tipp and Tipj^ act on the combined 
space of valley, sublattice, and spin degrees of freedom, 
and represent the normal and anomalous part of the BdG 
Hamiltonian, respectively. These blocks satisfy 



I i-pp / i-pp 



'Iph 



ph 



(Hermiticity) , 
(Fermi statistics). 



(9.2) 



To represent the single particle Hamiltonian /C, define 
the 256 16-dimensional Hermitian matrices 



X.. 



A"! 



P-3 



f^4 



(9.3) 



where /^i 2 3 4 = 0: 1? 2, 3. Here, we have introduced the 
four families p,, , s,, , a,, , and r,, of unit 2x2 and Pauli 
matrices that encode the particle- hole (Nambu), spin- 
1/2, valley, and sublattice degrees of freedom of graphene 
or the square lattice with 7r-flux phase, respectively. 

The Dirac kinetic energy /Cg of graphene or the square 
lattice wit h 7r-f lux phase that accounts for the BdG block 
structure (E).lc) is assigned the two 16x 16 Dirac matrices 



"1 = -'^oosi: "2 = -'^3032 > (9-4a) 
and is given by 

/Co := a • {-id). (9.4b) 

Similarly, by introducing the 16x16 Hermitian matrices 

/3 = ^3010' -^=^3033' 75 = ^3030' (9-4c) 

the counterpart to H in Eq. (2.1) is given by 

/C := /Co + /Cg^^g^ + /C.^g^i^,, (9.4d) 

where 



/Cg 



;alar 



a • (-a - 0575) , 



(9.4e) 



Given the Dirac kinetic term ICr,, we treat X,, „ ,, ,, 
as a perturbation, 

where m G M is constant in space and time. If X,, ,, ,, ,, 
anticommutes with the Dirac kinetic energy /Co, then it 
opens a gap in the massless Dirac spectrum of ICq. Wc 
shall call such a perturbation a mass in short. Each mass 



can be thought of as being induced by a breaking of a 
microscopic symmetry (see below). 

There are 64 = 4 x 16 mass matrices (i.e., X,, ,, ,, ,, 
that anticommutes with ICq). Of these 64 mass matrices, 
only 36 satisfy the condition 



-X 



(9.6) 



for particle-hole symmetry (PHS) and are thus compati- 
ble with the symmetry condition {p^ Sq(E)<Tq(^ ''0^)* — 



^'t on the Nambu spinors [i.e., compatible with Eq. (9.2)]. 
All mass matrices with PHS are enumerated in Table y. 

All 36 mass matrices from Table || can be classified in 
terms of the following (microscopic) 3 symmetry proper- 
ties, (i) A BdG Hamiltonian is time-reversal symmetry 
(TRS) when 



^0211 ^ ^0211 ~ ^- 



(9.7) 



(ii) A BdG Hamiltonian has SU(2) spin rotation symme- 
try (SRS) when 

[-'^31001^] = [-'^02001^] = [-'^3300)^] = 0. (9.8) 

(iii) A BdG Hamiltonian has sublattice symmetry (SLS) 
when 



^0033 ^ ^0033 



-/C. 



(9.9) 



For any lattice rcgularization of the BdG Hamilto- 



nian (9.5) supporting two sublattices A^ and Ag, as is 



the case for graphene or the square lattice with 7r-flux 
phase, the microscopic order parameter corresponding to 
a mass matrix satisfying the SLS (9.9) is a non- vanishing 
expectation value for a fermion bilinear with the two lat- 
tice fermions residing on the opposite ends of a bond con- 
necting a site belonging to sublattice A^ and another site 
belonging to sublattice Ag. We shall say that such a mass 
matrix is associated to a valence-bond solid (VBS) order 
parameter in analogy to the terminology used for quan- 
tum dimer models. A VBS order picks up a microscopic 
orientation that translates into a complex-valued order 
parameter in the continuum limit. Hence, we shall distin- 
guish between the real (ReVBS) and imaginary (ImVBS) 
parts of the VBS. Triplet superconductivity is also possi- 
ble on bonds connecting the two sublattices. The termi- 
nology TSC will then also be used. To distinguish TSC 
with or without TRS we shall reserve the prefixes Re and 
Im for real and imaginary parts. This is a different con- 
vention for the use of the prefixes Re and Im than for a 
VBS. 

Any mass matrix that does not satisfy the SLS 



9) 



corresponds to a microscopic order parameter for which 
the fermion bilinear has the two lattice fermions sitting 
on the same sublattice. Microscopic examples are charge- 
density waves (CDW), spin-density waves (SDW) such 
as Neel ordering, orbital currents leading to the quan- 
tum Hall effect (QHE), spin-orbit couplings leading to 
the quantum spin Hall effect (QSHE), singlet supercon- 
ductivity (SSC), or triplet superconductivity (TSC). 
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When SU(2) spin symmetry is broken by the order 
parameter, we add a subindex x, y, or z that specifies 
the relevant quantization axis to the name of the mass 
matrix. Moreover, TSC with SLS must be distinguished 
by the 2 possible bond orientations (the underlying two- 
dimensional lattice has 2 independent vectors connecting 
nearest- neighbor sites). These 2 orientations arc specified 
by the Pauli matrices used in the valley and sublatticc 
subspaces, i.e., by the 2 pairs of numbers 02 and 32. Sym- 
metry properties of all 36 PHS masses are summarized 
in Table ||. 

The set of all 36 PHS masses in Table g is invariant 
under an involutivc transformation defined by 



(9.10) 



and which we shall call C conjugation to distinguish it 
from the particle-hole transformation ( |9.6D . Here, = 
(s3±Sg)/2. For graphene or the square lattice with tt-Aux 
phase, this transformation corresponds to 



'•aT' 



'•fit' 



(9.11) 



where at » and bt „ creates an electron with spin s —1: 
, I on sublattice and sublattice Ag, respectively (see 



Appendix Under this transformation 



X 



(9.12) 



Hence, it leaves the massless Dirac kernel /Cg invariant. 

The organization of the mass matrices in Table ^ can 
be understood as follows. 

First, we preserve both SRS and charge conservation, 
i.e., we start with the 4 order parameters we have already 
encountered in the spinless case with charge conservation. 
There are two valence bond solids, ReVBS (Re A) and 
IniVBS (ImA). They have maximal symmetry a nd ar e 
invariant under the operation of C conjugation ( |9.12 ). 
The CDW order parameter (fi^) breaks the SLS. It is 
mapped into the Ncel spin-density wave with quantiza- 
tion axis z under the operation of C conjugation (9.12). 
The QHE order parameter (ry) breaks both the SLS and 
TRS symmetri es. It is invariant under the operation of 
C conjugation ( |9.12 ). 

Second, we break SRS with or without either TRS or 
SLS while always preserving charge conservation. The 
breaking of SRS is achieved by choosing a preferred quan- 
tization axis, say x, y, or z in SU(2) spin space. Breaking 
SRS while preserving SLS is achieved with spin-polarized 
valence-bond ordering in 6=3x2 different ways, which 
we abbreviate by ReVBS^, RcVBS^^, ReVBS^,'lmVBS^, 
ImVBSy, and ImVBS^ in Tabic |l| In doing so TRS is 
always broken. Breaking SRS and SLS while preserv- 
ing TRS is achieved through any of the 3 order parame- 
ters for the spin quantum Hall effect (QSHE) introduced 
by Kane and Mele in Rcf. E9|, which wc abbreviate by 



QSHE^, QSHEj^, and QSHE, in Table ||. Breaking SRS, 
SLS, and TRS is achieved through any one of 3 colin- 
ear magnetic order in the form of Neel order, which we 
abbreviate by Neel^, Neel^^, and Neel, in Table |l|. 

This brings the number of order parameters that con- 
serve the electronic charge to 16=4+6+3+3. There are 
thus 20=2+6+6+6 remaining order parameters that do 
not conserve the electronic charge. 

Third, superconducting order is achieved microscopi- 
cally by pairing two electrons sitting on different or iden- 
tical sublattices. In the former case, SLS is preserved. In 
the latter case, SLS is broken. Pairing of the 2 electronic 
spins takes place cither in a singlet or in a triplet channel. 
Antisymmetry under exchange of the two electrons mak- 
ing up a spin-singlet Cooper pair can only be achieved in 
an even angular momentum channel. On-site pairing is 
of course associated to vanishing angular momentum so 
that singlet superconductivity can only be realized when 
SLS is broken. This only leaves 2 possible singlet super- 
conducting order parameters that arc distinguished by 
whether they preserve or break TRS. They are denoted 
ReSSC and ImSSC, respectively. (Real and imaginary 
parts thus take a different meaning here as for ReVBS 
and ImVBS.) 

Fourth, a triplet superconducting order parameter, 
which we abbreviate by TSC in Table ||, is character- 
ized by a vector d in SU(2) spin space. This vector 
can point along any one of the three quantization axis 
x, y, and z in SU(2) spin space. Moreover, it can ei- 
ther preserve or break TRS for which cases we use the 
notations ReTSC and ImTSC, respectively, in Table ||. 
(Real and imaginary parts thus take a different mean- 
ing here as for ReVBS and ImVBS.) When SLS is pre- 
served by the superconducting order parameter, there 
are 12 = 2x2x3 independent order parameters, for a 
second factor of 2 besides the one for TRS arises since 
there are 2 directed nearest-neighbor lattice-bonds con- 
necting nearest-neighbor sites of the two-dimensional lat- 
tice. This is abbreviated in Table by using the index 
bond=02,32 in ReTSC,„„,„ ReTSC,,„,,^, ReTSC,,„,„ 



ImTSCu 



ImTSC^ 



and ImTSCu 



-'bondx' ^^^^^^'-^bondy '^^^'^ ^ ^'-"bondz • Finally, 

when SLS is broken by the superconducting order pa- 
rameter, there are 6 = 2x3 independent order parame- 
ters that we abbreviate by RcTSC^, ReTSC^^, ReTSC,, 
ImTSC^, ImTSCj^, and ImTSC, in Table 0. 

There are 12 = 4x3 order parameters th at are invari- 
ant under the operation of C conjugation (|9.12 ). They 
can be arranged in 4 groups of 3 each. Each group of 
3 obeys the same algebra. The 4 groups of 3 are: (i) 
ReVBS, ImVBS, QHE; (ii) ReVBS,, ImVBS,, QSHE,; 
(iii) ReTSCo2:r> ImTSC32^, ImTSC^; (iv) ReTSC, 
ImTSCgaj^, and ImTSC, 



02y: 



The operation of C conjugation ( 9.12 ) is a useful tool 
to identify the possibility of exotic topological effects. 

For example, we observe that the pair of SSC order 
parameters ReSSC and ImSSC, studied in Refs. |5^ and 
|5l| in the context of graphene, are conjugate by C to 
the Neel order parameters Neel^. and Neel^ , respectively. 
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Furthermore, Table |l| indicates that several triplets of 
masses that obeys the SU(2 ) alg ebra are related by the 
operation of C conjugation ( |9.12 ). They are 



(9.13) 



{ReVBS,ReVBS,CDW} 

^ {ReVBS, ImVBS, Neel J , 
{ReVBS^, RcVBSy, ReVBS J 

^ {ImTSC32„ ReTSC32„ ReVBS J , 
{imVBS.^, ImVBSy, ImVBS J 

^ {ReTSCo2„ImTSCo2.,ImVBSJ , 
{QSHE^, QSHEy, QSHEj 

^ {ImTSC^, ReTSC^, QSHEJ . 



Vortex-like defective textures in any of these mass dou- 
blets or meron-like defective textures in any of these mass 
triplets display fractionalization of some suitably defined 
quantum numbers. 

Finally, the topological property that a band insulator 
supporting the QSHE carries an odd number of Kramers 
doublets on its edges carries over to the C conjugate TSC. 
More precisely, the fact that the superconductors with 
the ImTSC^ and ReTSC^ order parameters are exam- 
ples of Zo topological triplet superconductors according 
to Refs. ^ |3[ IsJ, ^ and|5^ is here a mere consequence 
of their C conjugation with the QSHE^ and QSHE^^ order 
parameters, respectively. 



B. Classification of 5-tuplets of masses in 
graphene and vr-flux phases 

Mass matrices that commute pairwise generate com- 
peting local order parameters. Conversely, mass matri- 
ces that anticommute pairwise generate compatible local 
order parameters. 

All but one PHS masses anticommute with 16 out of 
the 36 PHS masses. The Haldane mass is unique in that 
it commutes with all PHS masses. 

There are 560 sets of three mutually anticommuting 
PHS masses. These triplets are generalizations of the 
triplet of compatible masses A = Re A -I- ilm A and /x^ . 
Integration over the Dirac fermions in the presence of any 
one of these mass triplets of mass m in competition with 
the Haldane mass ry induces an 0(3) NLSM in (2+1)- 
dimensional space and time with or without a Hopf term 
for m > \ri\ and \r]\ > m, respectively, as was derived in 
Ref. ^. 

There are 280 sets of four mutually anticommuting 
PHS masses and the maximum number of pairwise an- 
ticommuting PHS mass matrices is 5. Out of (^^^ = 
376992 possibilities, there are 56 distinct 5-tuplets of 
compatible PHS mass matrices. They are enumerated in 
Tables [II. (If PHS is not imposed, the maximum num- 
ber of pairwise anticommuting mass matrices in the 64 




FIG. 7: Two setups used to induce topological defects in an 
order parameter that support fractional quantum numbers. 



mass matrices is 7. There are 288 distinct 7-tuplets of 
compatible mass matrices.) 

In the background of each of these 5-tuplet, integra- 
tion over the Dirac fermions yields an 0(5) NLSM in 
(2+l)-dimcnsional space and time augmented by a Wess- 
Zumino-Witten (WZW) term as was derived in Refs. ^ 
and |5^. Defects-driven continuous phase transition be- 
tween phases of matter unrelated by symmetries (i.e.. 
Landau forbidden) become possible whenever the quan- 
tum numbers of the defective order parameters in a given 
5-tuplet jOre dual in the sense of BF Chern-Simons field 
theories.LJ We illustrate this idea with the following ex- 
amples. 



VBS-SSC-CDW 5-tuplet 



The 5-tuplet 



{ReVBS, ImVBS, ReSSC, ImSSC, CDW} (9.14) 

embeds the triplet made of the rGDW and the 2 VBS 
order parameters into a 5-tuplet .e3 Integration over the 
fermions yields an 0(5) NLSM augmented by a WZW 
term for the corresponding 5-tuplcts of bosonic fields 
111, 712, '^3 1 '^4j ^^'^ '^s obeying the constraint that 
they add in quadrature to unity. The 0(5) symmetry 
can be broken, either spontaneously or explicitly, down 
to the U(l)xU(l) subgroup corresponding to holding 



^CDW 



'51 



^BDW 



and A 



ssc 



fixed (except at the core of topological defects) through- 
out space and time. The corresponding Goldstone modes 
are the phases ^bd^ and OggQ- They become charge 
2 Higgs fields if the U(l)xU(l) global symmetry they 
generate is gauged through the introduction of the ax- 
ial gauge fields a'^ 



ygg and the electro-magnetic gauge 
Their dynamics i s gov erned 
by the Andcrson-Higgs-Chcrn-Simons theory (4.13) with 



fields OggQ, respectively. 



the identifications 



^BDW' 



as 



^VBS' 



and 



^ssc ~ "^^^ VBS phase is destroyed when 

the vortices carried by the conserved topological current 
j^^i^ = e^''''9^ap6'vBs/(27r) deconfine. The SSC phase 
is destroyed when the vortices carried by the conserved 
topological current jgs*c = ^^"^(^u'^p^ssd (2'"') deconfine. 
Because of the BF term in the effective action, the quasi- 
particlcs supported by JvBS ^^^"^ carry a fraction of the 
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gauge charge of the gauge fields a^g^ , while the quasipar- 
ticles supported by igg*^ also carry a fraction of the gauge 
charge of the gauge fields ayBS- Furthermore, both types 
of quasiparticles are bosons (there is no TRS-breaking 
Haldane mass) . From these two facts follows that decon- 
finement of one type of quasiparticles implies confinement 
of the second type of quasiparticles, i.e., a direct transi- 
tion between the VBS and SSC phases. 

An experimental setup to detect exotic quantum nuiii=. 
bers related to the 5-tuplet (9.14) is given in Fig. |^(a).L2l 
We assume that graphene sits on top of a type-II s-wave 
SC substrate. By the proximity effect, graphene develops 
a SSC order. The SSC order can coexist with the C DW 
and VBS orders in graphene according to Eq. (9.14). An 



applied magnetic field perpendicular to graphene creates 
an Abrikosov lattice of vortices in the substrate and, by 
the proximity effect, in graphene. The magnetic fiux 
tube s thre ading graphene pin axial charges according to 
Eq. ( [4.13[ ). (See also Refs. |l| and ||.) Increasing the 
magnetic field so as to destroy SSC deconfines the axial 
charges, i.e., stabilizes the VBS. Conversely, destroying 
the VBS by the deconfinement of VBS vortices also de- 
confines the electric charges, i.e., stabilizes the SSC. 



An experimental setup to detect exotic quantum num- 
bers related to the 5-tuplet ( 9.16 ) is given in Fig. 0(b). 
We bring in contact a (3D) bulk type-II SSC with a mate- 
rial displaying the QSHE. Instead of graphene for which 
the spin-orbit coupling is very small, HgTej/fHfrrCd)Te 
semiconductor quantum wells are suitable.OOO Any 
SSC vortex in the substrate induces by proximity effect 
an "5^ spin charge" in the device supporting the QSHE, 
while any "spin flux" in the device supporting the 
QSHE induces an electric charge. 



4. XY-Neel-TSC-QSHE 5-tuplet 



The oper ation of C conjueatiop ( 9.12 ) on the 5- 
tuplet (|9!l6|) yields the 5-tupletEl'L3 



{Neel^, Neelj^, ImTSC^, ReTSC^, QSHE J. (9.17a) 

By rotating SU(2) spin quantization axis, i.e., by cyclic 
permutation of the indices x, y, and z, we also get the 
5-tuplets 

{NeeL, Neel^, ImTSC^, RcTSC^, QSHE^} (9.17b) 



VBS-Neel 5-tuplet 



and 



The oper ation of C conjugation ( 9.12| ) on the 5- 
tuplet ( p.l^ yields the 5-tuplet 



{ReVBS, ImVBS, Neel^, NeeL, Neel J. 



(9.15) 



The triplet of Neel order parameters is here embed- 
ded into a 5-tupl et .. by adding the doublet of VBS or- 
der parameters.E3E§E2l This 5-tuplet has been discussed 
in the context of deconfined quantum criticality of two- 

rlimpnsinrial .S* 1 /? qnantnm antifeiTOmagnetic Spin 

models-EyaOElHO The 5-tuplet (|I|) is the only 
5-tuplet supporting the full SU(2) symme try of t he Neel 
vector. The symmetry analysis of Sec. IX B 1 follows 
with the identifications O^rvta — > ^vbS' ^ssc ^ ^Ncoi ' 

^VBS ~^ JVBS' ^^'^ 



"VBS 
■ vrt/i 

Jssc ' 



^VBS' "^SSC 



^VBS 



■ vrtfi 



SSC-QSHE 5-tuplet 



{Neel^, Neel^, ImTSC^, RcTSC^^, QSHEJ. (9.17c) 

These 5-tuplets describe SLS- and SRS-breaking order 
parameters consisting of an easy plane antiferromagnetic 
order parameter coexisting with the QSHE and T SC or- 
dcr parameters. The symmetry analysis of Sec. llXBl 
follows with the identifications 



"TSC^' '^VBS 
J -vrtu 

and jgs^ 



^NccU 



^SSC 



^VBS 



^TSC^ 



-vrt/^ 
■?VBS 



-vrt/^ 



-vrt^ 
■?TSC^ 



say. 



An experimental setup to dete ct ex otic quantum num- 
bers related to the 5-tuplet ( 9.17| ) is also given in 
Fig. 0(b). Any defect in the bulk XY antiferromagnct, 
i.e., a magnetic vortex, induces a localized midgap state 
that carries a fraction of the electric charge carried by 
the phase of the TSC in the band insulator supporting 
the QSHE. Any TSC vortex induces an "5'^ spin charge" 
in the device supporting the QSHE. [A related fractional 
(electrical) charge is discussed at the helical edges of the 
QSHEO] 



The 5-tuplet 



X. 



DISCUSSION 



{ReSSC, ImSSC, QSHE^, QSHE^, QSHE J (9.16) 

embeds the triplet of QSHE order parameters into a 5- 
tupl et || by adding the two possible SSC order p arame- 
ters.BE^O The symmetry analysis of Sec. IX B 1 follows 
with the identifications 0, 



^VBS 



a. 



QSHE^ 



-vrt/^ 

■?VBS ~^ 



Jssc'' ^'^'^ 



■ vrt^ 

Jssc 



-vrt/^ 
■?QSHE, 



Motivated by the interplay between charge-density 
(/ig), bond-density (A = |A|e~'^), and integer quantum 
Hall (77) instabilities in graphene-like two-dimensional 
electronic systems, we have computed the fractional 
charge and fractional statistics of both screened and un- 
screened quasiparticles. 

At the microscopic level, screened quasiparticles are 
here the linear superpositions of two bond-density waves 



25 



(A and Og), each of which carry a point defect. Un- 
screened quasiparticles are defects in one type (A) of 
bond-density wave. 

In the long-wave-length and low-energy limit and af- 
ter integrating out the fermions, the quantum dynamics 
of scree ned quasiparticles is controlled by the effective 
theory (4.13) of the Anderson-Higgs-Chern-Simons type 
involving three fields. There arc two U(l) gauge fields 
and one phase field. 

The first gauge field is responsible for the conser- 
vation of the total fermion number. The second gauge 
field is responsible for the conservation of a relative 
fermion number, i.e., the difference in the fermion num- 
ber located at the two valleys of graphene say, and is thus 
called an axial gauge field. The phase field 6 ~ — argA 
originates microscopically from the fact that bond distor- 
tions include atomic displacements away from the crys- 
talline order that are parametrized by continuous angular 
degrees of freedom. 

Screened quasiparticles are not yet explicitly mani- 
fest in the field theory ( 4.15 ). They appear as point 
particles with the conserved topological current j^^.^ = 
{2TT)~^€^^'^^djjd^6 that carries no axial gauge charge, once 
a duality transformation has been p erform ed. The La- 
grangian dual to the Lagrangian ( [1.13 ) can be pre- 
sented as aJZIhern-Simons theory for 4 gauge fields whose 
K matrixllj is 4-dimensional and couple through a 4- 
dimensional charge vector to the vortex cwrent. Because 
the K matrix has a vanishing eigenvaluc|[3 this dual the- 
ory is not a topoLaaical theory, say such as a BF Chern- 
Simons theory.ac3 The vanishing eigenvalue of the K 
matrix signals the existence of low-energy excitations, 
the screened quasiparticles. Their fractional charges Q 
and statistical angle Q can then be calculated and are 
presented in the phase diagram of Fig. 

When the U(l)xU(l) local gauge symmetry holds, 
i.e., for screened quasiparticles that represent vortices in 
the phase field 9 whose axial charges are dynamically 
screened by axial gauge half fluxes in 05^, the fractional 
charge Q and the fractional statistical angle 8 in Fig. |l| 
are complementary. One is non-vanishing if and only if 
the other vanishes. Moreover, Q and Q are universal in 
the fully gaped phases for which they arc non-vanishing 
and given by a rational number in some units. 

When the U(l)xU(l) local gauge symmetry is bro- 
ken, i.e., for unscreened quasiparticles that represent vor- 
tices in the phase field 9 without the attachment of ax- 
ial gauge half fluxes, the fractional statistical angle 8 is 
non- vanishing everywhere in Fig. |l| with a discontinuous 
jump at TO = fJ.s + I Ap = It^I and a non-universal de- 
pendence on the ratios /^s/to and rj/m. The fractional 
charge Q is only non- vanishing when |?/| < to where it is 
also non-universal. 

Comparing the values of 8 in Fig. |l| calculated from 
field theory with a numerical evaluation of 8 for an un- 
derlying microscopic (lattice) model is difficult for two 
reasons. 

Defects in the phase 9 have a characteristic size of the 



order of 1/to for lattice models, i.e., they bind a fcrmionic 
charge through midgap states. The profile of defects in 
the axial gauge fields Og^ is power law, i.e., they bind 
a fcrmionic charge through threshold continuum states. 
Thus, the linear extend of any lattice model must be 
much larger than 1 /m for any reliable numerical calcula- 
tion of 8. On the one hand, if we impose the U(l)xU(l) 
local gauge invariancc at the lattice level, the system sizes 
accessible to a numerical computation of 8 are, at best, 
of the order 1/m, i.e., too small for a comparison with 
field theory. On the other hand, if the U(l)xU(l) local 
gauge invariance does not hold at the lattice level, say af- 
ter performing a mean-field approximation for which the 
accessible system sizes are sufficient to measure Q with 
the help of a static probe such as the spectral asymmetry, 
then the values of Q and 8 arc not universal anymore. 
To put it differently, the values of Q and 8 measured 
dynamically depend sensitively on the dynamical rules 
used. But these dynamical rules are model dependent 
when they are not fixed by imposing the local axial gauge 
symmetry. 

The fractional charge Q or the statistical angle 8 in 
the phase diagram of Fig. |l| disagree with the results of 
Rcfs. |l| 0, and |l|. 

Although the charge assignment in Ref. [l^ agrees with 
that in Fig. |^ the statistical angle ia ascribed the value 
8 = sgn(?/)7r/4 whenever 77 7^ 0.L3 However, the sta- 
tistical angle 8 is non- vanishing if and only if the Hopf 
term is present in the 0(3) non-linear-sigma model de- 
rived in Ref. |l^, i.e., if and only if I77I > to, in which 
case full agreement with the charge and statistical angle 
assignments of Fig. ^ is recovered. 

Seradjch and Franz in Ref. |l^ have computed the 
fractional charge Q and fractional statistics 8 of dy- 
namical defects in 8 and for the field theory (3.1) 
when = ?7 = 0. Their analysis has been repeated 
by Milovanovic in Ref. |l8|. They found the assignments 
Q = ±1/2 and 8 = ±7r. Their scmion statistics contra- 
dicts our result 8 = in Fig. |l]. This discrepancy can be 
traced to the fact that Seradjch and Franz used a singular 
chiral U(l) gauge transformation with the Pauli-Villars 
regu lariza tion to derive an effective action different than 
Eq. (4.13). As we show in Appendix ^ the effective action 
used by Seradjch and Franz, when suitably generalized 
to the case fi^ ^ rj = 0, fails to reproduce the fractional 
charge (5.7) of quasiparticles in the presence of a flux in 
gauge field. Explicitly, it follows from Eq. 9 of their 
paper Ref. |l^ that the fractional charge in the case when 
the mass vortex is accompanied by an axial half-fiux, en- 
forcing the screening condition — ^d^9 = 0, is Q = ! 
However, the fractional charge (5 = 1/2 [see Eq. (^?7|)] of 
screened quasiparticles is a result established from direct 
(static) numerical computation of Q on a suitable lattice 
regularization of the field theory (3.1). 

The charge-density {fig), bond-density (A), and inte- 
ger quantum Hall (77) instabilities are the only instabil- 
ities compatible with the electron-number conservation 
and SU(2) spin-rotation symmetry (these are, naturally. 
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also the only four possible instabilities for the spinless 
case). However, there can also be superconducting insta- 
bilities or. if the electron spin is accounted for, magnetic 
instabilities. We have performed a systematic classifica- 
tion of all instabilities for the f6 dimensional free Dirac 
Hamiltonian induced by local order parameters that re- 
spect the Bogoliubov-de-Gcnnes particle-hole symmetry. 
We have found that the order parameter for the integer 
quantum Hall effect (Haldane mass 77) is unique, for it 
competes with all other instabilities. We have also found 
that the largest number of coexisting order parameters 
is 5 and enumerated all the corresponding 5-tuplets of 
masses. Each of these 5-tuplet can be thought of as a gen- 
eralization of the 3-tuplet (/Xg,Re A,ImA) that supports 
quasiparticles with fractional quantum numbers. These 
5-tuplcts provide a rich playground for Landau-forbidden 
continuous phase transitions. Any U(l) order parameter 
in a 5-tuplet can be assigned a conserved charge and sup- 
ports topological defects in the form of vortices. A pair 
of U(l) order parameters in a 5-tuplet is said to be dual 
if the vortices of one order parameter binds the charge 
of the other order parameter and vice versa. A contin- 
uous phase transition can then connect directly the two 
dual U(l) ordered phases through a confining-dcconfining 
transition of their vortices. 
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APPENDIX A: CALCULATIONS OF THE 
COEFFICIENTS Cj?\ C^,V, C^W C[l\ AND 0,^3' 

Let be a 4-dimensional representation of an element 
of the Lie Algebra u(2) generated by the unit 4x4 matrix 
and the 4x4 matrices T,^ = (S^, Sj, S3) whereby 



Pa, ^b] = ^fabc^c: {^a: ^b} = 2(5, 



ab' 



(Al) 



In this appendix, we are going to integrate the Grass- 
mann fields x ^^'^ X in the partition function 



S := I dx^dx^dx^C 



(A2) 



Here, the Feynman slash notation 



(A3a) 



is used when contracting 3- vectors with the 4x4 matrices 
F^ = (Fg,F]^,F2) that realize a 4-dimensional represen- 
tation of the algebra 

{F^,Fj = 2g^„ .9^, =diag(l,-l,-l), (A3b) 

while they commute with Y,^ = (Sj^, Sj, Sg), 



(A3c) 



We shall work in momentum space. To this end, we 
introduce the Fourier transforms 



X{x) = 
X{x) = 



Ak-x 



Xik), 



(A4a) 



whereby the notations 

dk°dk^dP 
(27r)3 ' 



(A4b) 



will be used. The action and Lagrangian in Eq. ( A2) are 
represented in momentum space by 



S = 



The free propagator, here defined by 
1 



1 - ^2)] xih)- 

(A5) 



Goik) := - 



^ — 77 — mEg 



^2 _ _ ^2 _ 2r]mT,^ 

{ft — rj — mTi^) (/c^ — jf' — ni^ + 2rymSg) 

(A6a) 

can be decomposed into the sum of the unit 4x4 matrix 
weighted by the factor P{k) and the 4x4 matrix Sg 
weighted by the factor Q(fc); 

Go(fc)=F(fc) + g(fc) Sg, 

^ ik"^ - if - w?) -riik"^ - if + m^) 



P{k) 
Q{k) 



(fc2 _ _ rnf^ (k^ - (77 + m)^) 



(A6b) 
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The induced effective action for the background If} is 
defined by 



(A7) 



?i=i 



where it is understood that the Grassmann integration is 
performed in a w ay th at preserves the local U(l) xU(l ) 
gauge symmetry ( 4.3b| ). The effective action in Eq. ( 4.1(]| ) 
with the coefficients from Table | follow s by combining 
the local U(l) xU(l) gauge symmetry (Ob) with the 
loop expansion (A7) up to the order n = 2, 



tr [Go(fc)^(0)] , 

tr [Go(fc)^(q)Go(fc-g)^(-g)], 



(2)r 



(A8) 



One verifies by explicit calculation that 



(A9) 



To proceed wit h the evalua tion o f S^^^ , we note that 
the algebra Q, ( |A3b| ), and ( |A3^ ) can always be realized 
with the choice 



7» 



1. 



1, 



(AlO) 



where 7^ = (7o>7i,72) and = (0-^,(72,(73) realize two- 
dimensional representations of the algebra ( A3b ) and 
(Al), respectively. With this choice, it is obvious that 
a single trace over the 4x4 matrices spanned by the unit 
4 X 4 matrix, = (ro,ri,r2), and = {Y.^,T.^,T.^) 
factorizes into the product over two traces; one trace over 
the 2x2 matrices spanned by the unit matrix 1^ and 
7ai ~ (7oi7i'72) ai^d one trace over the 2x2 matrices 
spanned by the unit matrix 1^ and = (ctj^, ctj, (Tg). It 
then follows that 



- I [PQiq)rtTABAqW3B,{~q)] 



where [H K (q)]'^'^ , H and K being P or Q, is defined as 



[HKiq)]" 



tr[Hik)YKik-q)-f^], (Allb) 



and with the understanding that 



B.. 



(Allc) 



If the integrals ( |Allb| ) are regularized so as to preserve 
the Lorentz covariance, then they must be of the form 



[HKiq)r = g-"^ [HK] 



(0) 



(A12) 



to linear order in q. Fu rtherm ore, imposing a regular- 
ization of the integrals (Allb ) th at preserves the local 
U(l)xU(l) gauge symmetry ( 1.3b ) demands that the co- 
efficients 



•-^00 



•-^33 



C, 



[PPf^ + [QQf^ 

■.=2{[PQr+[Qpr 

This gives the effective Lagrangian 



(A13) 



C. 



(2) 
eft 



+ Cl})6'''"' {bld^bl + bld^bl + ...) 



(A14) 



C^^e'^P'^bldJil 



with the coefficients 



(0) 



1(0) 



•-^00 



(1) 

33 



[PPf' + [QQ] 



(1) 



C« := [PPr - [QQ] 



(1) 



1(1) 



(A15) 



■.^2[[PQr + [QP] 



(1) 



[QQ{q)Y''trA^Mq)^zB,{-q)] 
The coefficients ( |Al5| ) are evaluated by pcrforn^^i^-|^a^Wick rotation to the Euclidean metric with the rules 



Under these rules 



7", 



1 



ui/A 



a = 1,2,3. 



(A16) 



(A17) 
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while the scalar functions P{k) and Q(k) in the propagator ( A6b| ) take the form 

P{k) = [i?(fc)^ + Sik)] ® 1^, Q{k) = [Tik)^ + C/(fc)] ® a„ 

with m_|_ := m±r] and 

z ( fc^ + 



(Al8a) 



Rik) 
T{k) = 



2 m^+^Ti'^ 



(fc2 +m2_) (fc2 +to2_) 

— 27]mi 



t/(fc) 



— 77 (fc2 — ?7!,-|_?7!,_) 

— m (fc2 + m-|_?7i_)) 
(fc2 + ?71^) (fc2 + mi) ' 



(AlSb) 



(fc2+m2_) (fc2+771^)' 

Their small q expansion are 

R{k -q)= R.{k) + (k ■ q)R^^'> (k), S{k -q)^ S{k) + {k ■ q)S^^^ (k), 



where 



i?(i)(fc) 
S'^'\k) 
T'^^\k) 
U^^\k) 



T{k -q)= T{k) + [k ■ q)T^^\k), U{k -q) = U{k) + [k ■ q)U^^\k), 



(A19a) 



2i 



(fc2- 


1- m\Y (fc^ - 
-2r; 


h m2^)^ 


(A:2^ 


- m\Y {k^ 4 


- m2_)^ 


(fc2H 


-2m 


-m2_)' 


(F- 


F m^)^ (/c2 - 





l^fc'' + {m\ + m^) fc^ - m\Tn^_ + 
[k^ - 2k^m^ 



^mj_m_ — mim?_ — m,m_ imi_ + rn?_ 



(A19b) 



[2fc2 + (m^ + m,2,)] , 

[fc'^ + 2A:^m_|_m_ — m.^m^ + m_^m_ (m^ + m.^) 
I 



At last, the coefficients (A15) follow from 



C 



(0) 

11 



cg^^cg^=-^{[PPr + [QQ] 



^00 



l(i) 



l(i) 



(1) 

03 



27. (\PQf^ + [QPf^ 



(A20a) 



with 



[QPf^ 
[QPf^ 



in 
(2^ 

Btt 
(2^ 



(A20d) 



fpp](0) 

[PP]{1) 



47r 
(2^ 

(2^ 



47r 
(2^ 

Stt 
(2^ 



"^^RSW ^ ^SR + g-^Sfld) ; ' 



(A20b) 



47r 
(2^ 

Stt 
(27r)3 

Here, the integrals 
/ 



(A20c) 



CXD 



f ~ 1 



XY = J dkk^X{k)Y(k), 


ixY = j dkk^X{k)Y{k), 




(A20f) 



(A20c) 



with X{k) and Y{k) denoting i?(fc), i?(i'(fc), S{k), 
5(i)(fc), T(fc), T(i)(fc), f/(fc), or U^^\k) need to be regu- 
larized in a way that preserves the local U(f )xU(l) gauge 
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symmetry ( |4.3b| ) and the Lorcntz covariance. The brute 
force method consisting in imposing the ultra-violet cut- 
off A in the integrals (A20f) and ignoring all the terms 
linear in A delivers the coefficients from Table |[ 



APPENDIX B: DUALITY AND STATISTICS IN 
THE QUANTUM XY MODEL WITH A 
CHERN-SIMONS TERM 



The presentation of the effective action (4.13) with the 
help of Table | is not optimal for the purpose of extract- 
ing the statistical angle Q acquired by the pairwise ex- 
change of unit vortices from Sec. ^ Needed is a con- 
served vortex current that accounts for the local vortex 
density and the local vortex current generated by the 
physical process involving the exchange of two vortices. 
This vortex current can be non-vanishing anywhere in 
the phase diagram in Fig. ^. Thus, an optimal presen- 
tation of the effective action ( [4.13| ) should include this 
vortex current. This can be achieved by taking advan- 
tage of the duality between the quantum XY model and 
compact quantum [GlcGlpttdynamics in (2-|-l)-dimensional 
space and timeJI3'L3li3'EJH which we now briefly adapt 
for our purpose. 



We wish to constrain all configurations 9 appearing in 
the partition function Z^y by the condition (B4), i.e., 
we wish to restrict 6 to any configuration such that it 
supports the conserved current j^j.^ . We call such config- 
urations vortex configurations. 



The condition of current conservation (B4) can be en- 
forced by the three Lagrange multipliers ^ ~ 0, f,2. 
If so, the following partition function restricted to vortex 
configurations follows. 



(B5a) 



with 



-XY vrt 



K 



+ — e^^^d^d.dx + nd^ (B5b) 
47r 

I 



— 7^ dx9 < ^fj. . 

ZTT 



'^x" + Cfj, .7 vrt 

Here, we have introduced the field 



e^^f^d c 



(B6a) 



1. Duality 

We begin by defining the partition function for the 
quantum XY model in (2-|-l)-dimensional space and time 
with an additional Chern-Simons term. 



7CS 



J^xyW (5,^) (d^O), 



(Bl) 



The Chern-Simons action C(jg can be rewritten using an 
auxiliary vector gauge field d^, 

Z^^ J V[e] V[d^] e»/rf'^^(£x^[e]+c.[e,dM]) (B2a) 

with 

C,[9, d^] := ^ e^'^^d^d^d^ + ^ d^e^''^ 8,8x9. (B2b) 

Observe that the Chern-Simons gauge field d^ couples to 
the current 



= ^e^-^^8,dx9. 



Jvrt 



27r 



This current is necessarily conserved 



= 9^ J vrt • 



(B3) 



(B4) 



whose dual field is given by 



8 c 



(B6b) 



and we dropped total derivatives after performing partial 
integrations. 

The equation of motion for 9 gives the condition 



8^'9 



1 



2ttK 



7*F, 



(B7) 



from which we recover the inhomogeneous Maxwell equa- 
tion 



7^ - 

■Jvrt ^ 



1 



(B8) 



After integration over 9, the partition function (B5) 
becomes 



^XYvit 



(B9a) 



with 



C'^^lc d 1 



47r 



flJvTt 

fiJ vrt ' 



(B9b) 



The dynamical gauge fields and d^ that couple to the 
vortex current j^rt have a Maxwell for and Chern- 
Simons for d^ kinetic energy. They endow the quantum 
theory with an explicit U(l) xU(l) local gauge symmetry. 
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Exchange statistics 



We turn our attention to the computation of the ex- 
change statistics of vortices with current j'^j.^ interacting 



through the Chern-Simons action 



47r 



vrt J 



(BIO) 



The relationship between the current and the field that 
results from the equations of motion is 



Ztt 



Hence, the vorticity 



(Bll) 



(B12) 



supported by the vortex current j^j.^. is related to the 
circulation from the gauge potential through 

ng^-^J d^r {d.d^ - d^d,) ^-1-j^dl-d. (B13) 

Consider now winding two vortices, with vorticities n-^ 
and rij around each other. Without loss of generality, 
suppose that we hold vortex 1 at the location Xi fixed 
and move vortex 2 along any closed trajectory X2{t) that 

encircles once x^. On the one hand, the field dl!:\x2) at 
the location of vortex 2 that is indu ced by vortex 1 must 
then satisfy, according to Eq. (B13), 



dx2 ■ d^^\x2) = -27rni. (B14) 



On the other hand, the vector current resulting from 
moving vortex 2 around vortex 1 is 



Jvrt (''I ^) 



n^^Six-x^it)). (B15) 



Finally, the Berry phase acquired by winding vortex 2 
around vortex 1 is 



26 = — I 47m, 



dtS^^-'^ 

dt 



— ( 47m. 
47r 



2 ^ dX2 ■ S^' 



(B16) 



47r ^ 



We conclude that the statistical phase 8, which is one- 
half of the Berry phase, is given by 



e 

TT 



rii 712 



(B17) 



In particular, for a positive unit vortex winding around a 
negative unit vortex (anti- vortex), we find the statistical 
phase 



9 

— = K,. 
TT 



(B18) 



APPENDIX C: BERRY PHASE IN THE 
SINGLE-PARTICLE APPROXIMATION 

We are going to describe how Aharonov-Bohm phases 
7 or, more generally. Berry phases O accumulated un- 
der the pairwise exchanges of quasi-particles can be com- 
puted for non-interacting models of fermions defined on 
lattices. 

We first discuss Berry phases for lattice models of non- 
interacting fermions in all generality. We then specialize 
to the case of the tt flux phase for which we define vor- 
tices, axial gauge half fluxes, etc. 



1. Berry phase on the lattice 

Assume that we are given a lattice model, whose sites 
r and internal degrees of freedom are collectively denoted 
by the latin index m, that describes the quantum dynam- 
ics of non-interacting fermions. In second quantization, 
if the creation and annihilation c„ obey the usual 
fermion algebra 

{^m''^n} ^m,m {'^In^^li} {'^ini^nl (Cla) 

then we take our non-interacting Hamiltonian to be 

•■= -"^t^nnClnCn (Clb) 

where the matrix t with the matrix elements t„„ is Her- 
mitian. 



f =f 

mn nm' 



(Clc) 



We shall call the matrix t the background. Its uni- 
form diagonal matrix elements (the chemical potential) 
fixes the average number of fermions. We shall assume 
that some choices for the matrix t can be associated with 
point-like defects. These point-like defects can thus be 
labeled by their positions r^^ r2, ... on the lattice with 
their corresponding backgrounds t^^ ,^^ . For a given 
filling fraction, the many-body ground state in the back- 
ground ^ of point-like defects is the Fermi sea 



n 



(C2) 



Here, the prime over the product means that only the 
lowest single-particle energy eigcnstates arc to be filled up 
to the given filling fraction out of the state |0) annihilated 
by the c„. 

Imagine that we move the fc-th point-like defect along 
a closed path Vf, of counterclockwise orientation while 
holding all other point-like defects fixed. We then dis- 
cretize the path, thereby defining backgrounds , 
n = l N. 
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The gauge invariant phase 7^ is defined by 



N 



-^arg 



ri=l 



(n+1) 



(C3) 



If we do this exercise for two cases, one when the path 
Vj^ encircles another defect I and another one when the 
defect I hes outside the path V^, we obtain the statistical 
phase acquired by the counterclockwise exchange of 
point-like defects k and / from 



1 



T I inside T / outside 'P, 



(C4) 



This phase does not depend on the presence of other 
static point-defect inside the path V/^, for their con- 
tributions to 7 1 jjjgjjjQ -p cancel their contributions to 

^ I outside ' 

The overlaps 



(") 



An+1) 



(C5) 



from Eq. (C3) can be presented as the determinants for 
the products between two matrices built out of the eigen- 
vectors oft'"' and as we now show. For any back- 
ground t, define the unitary transformation U by 



UtU'^ = diag (e„ 



(C6) 



i.e.; U is the matrix of eigenvectors with energies of 
the single-particle Hermitian matrix t. For the two back- 
grounds entering the overlap ( |C5|) , these unitary trans- 
formations are denoted by and C/„_|_i, respectively. 
One then verifies that 



det {UiU„, 



(C7) 



Evaluation of the phases (C3) or ( |C4|) requires N diag- 
onalizations and the multiplication of N determinants, 
a computing exercise that scales as a power law in the 
number of sites in the lattice. 



The square lattice with a flux of tt per plaquettc (the tt 
flux phase in short) is the non-interacting tight-binding 
Hamiltonian 



- E E (^l!^+.,^i+.,«r. + H.c.) (C9a) 

with the (gauge dependent) choice of the tunneling am- 
plitudes 



Jtt/2 



t = it, 



^r,r+s., r.r-\-s. 



(C9b) 



Time-reversal symmetry is the property that H* is lo- 
cally gauge equivalent to H^. Sublattice symmetry is the 
property that —f — ff^ under the local gauge trans- 
formation 



r+s. 



(CIO) 



At half-filling, the Fermi surface collapses to two non- 
equivalent Fermi points due to the breaking of translation 
invariance, for the unit cell is now the unit cell of the 
sublattice with two atoms per unit cell. At half- 
filling, there are four non-equivalent ways to open a gap. 

There is the charge-density wave instability through 
the perturbation 



H.. 



tfj-s E v^^- 



(Cll) 



that breaks the sublattice symmetry of 
Hamiltonian (|C9|) under the local gauge transforma- 



tion (CIO) but preserves time- reversal symmetry, 



There is the bond-density wave instability through the 
perturbation 



2. Lattice defects for the tt flux phase 

Consider a square Bravais lattice A that is spanned 
by the orthogonal basis of vectors Si and of length 
a, the lattice spacing. We shall also define S3 = —s^ 
and S4 = — The square lattice is the union of two 
interpenetrating square lattices A^ and Ag with lattice 
spacing a. Any site rg S Ag can be decomposed in a 
unique way according to = + Sj^ with G A^. 

Because of the bipartite nature of the square lattice, we 
introduce fermionic annihilation operators denoted 

'a 

and h^^ and their adjoints for any site and of the 
sublattice A^ and Ag, respectively. These operators obey 
the usual fermionic algebra with the only non- vanishing 
anticommutators 



{a.^,StJ = l, {6^^,StJ = l. 



(C8) 



:= - E E [^^r,r+s^ bi+s^ K + H.C.) (C12a) 



with the tunneling amplitudes 
t 



],3J + (^J,2 + Sja)] 



X (A e+'i^e+"^-'' + c.c) 
where the wave vector 



K 



TT I 

a [0 



connects the two Fermi points 



^ 2a I 1 



(C12b) 



(C12c) 



(C12d) 
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[Here, G ■ r = + m2 if r = [m^s^ + 171282) with 
and ttij integers.] It preserves the sublattice and 
time-rcvcrsal symmetries of . Notice that St^ 

are purely imaginary (real) when j = 1,3 {j = 2,4). 
When the complex-valued order parameter A is turned 
into a space-dependent order parameter A = ^^{r) e'-'''''' 
trough an amph tude Ag(r) and phase /(r) modulation, 
then Eq. (|C12h| ) turns into 



iAo(r) 



> {5j.l+S,,s) + {S,,2 + S,,4)] 



X cos (/(r) + -j + G-r 



If the bond-density wave supports the unit vortex A(r) = 
Ao(r)e±*^ at the origin of the lattice, whereby we have 
introduced the polar coordinates r • Si/a = rcosO and 
r • 82/ a = rsmO, then 



St 



^[z(5,,+^,,3) + (^,,. + ^,,4)] ^^^^^ 

^ (C14) 



r,r+s . 



X cos [±9 + — j + G ■ r 

The case of an arbitrary distribution of vortices of integer 
charges n^'^' at the sites r*^*^^ follows with the identifica- 
tions 



/(r) = ^ n^'^^ arctan ■ 



„(fc) 



2 ' 2 



(C15) 



There is the time-reversal and sublattice symmetry- 
breaking bond-density wave 



reA. j=± 



E E(<.^t+a,^.+H.C. 
r6A„ j=± 



(C16a) 



where = Sj^ ± S2 



^2,- ~ ^2, 



t. (C16b) 



The lattice origin of the axial gauge field can also be 
identified with a staggered modulation of the nearest- 
neighbor hopping through the perturbation 



^5 := - E E [^tri+s^ hl+s^ + H.c.) (C17a 



with 



(5t 



(5) 



-5t 



(5) 



Motivated by the axial gauge flux 



O'Zii'^) = -na{r) e 



(C17b) 



(C18) 



where a(r) is any function that vanishes no slower than r 
at the origin and saturates to 1/2 at infinity that screens 
a charge n vortex in the continuum limit, we identify the 
lattice axial gauge flux that screens a charge n vortex 
located at the origin with 



6t 



(5) 



-6t 



(5) 



t—^ sin 61, 

r 



St^^l^^ =-^^1^1^- =-it^^^cos( 



(C19a) 



(C13) where we choose to regularize the vortex with 



air) = — tanh — . 
^ ' 2 C 



(C19b) 



Here, ^ is a characteristic length scale that determines the 
core radius of the axial gauge flux. The function a{r) reg- 
ularizes the singularity of 1/r at the origin. The case of a 
distribution of axial gauge fluxes located at r^''^ follows 
with the substitutions n — s- n'*^^ for the integer vortex 
charges, r — > r — r^*"'' for the positions of the axial gauge 
fluxes, and a linear superposition of the corresponding 
tunneling amplitudes. 

Finally, a uniform magnetic flux or a magnetic flux 
localized to one plaquette of the square lattice follows 
from the Peierls substitution 



r+s- 



(C20) 



in Eq. (C9b) with any suitable choice for the phases 



APPENDIX D: CHIRAL SINGULAR GAUGE 
TRANSFORMATIONS 



The effective theory (4.13) is one of the main results 
of this paper. From it follows the charge and statistics of 
quasiparticles. This effective theory was de rived b y com- 
bining symmetry arguments to reach Eq. (4.10b) and a 
direct computation to fix the coefficients that symmetry 
does not determine. Computation of these coefficients 
can be achieved in many independent ways. For exam- 
ple, the computation of the coefficient Cq^^ is fixed by 

obtaining the charge of quasiparticles. Hence, Cgg^ can 
be deduced from Rcfs. |l^. and |l^ for some range of 
the parameters /ig and rj, or. more directly, from numer- 
ics. The key step to derive the effective theo ry (4.13) 
was t he U (2) pure gauge transformation (4.1) [see also 
Eq. (3.1C)]. In this Appendix, we compare these dif- 



ferent ways of deriving effective actions for computing 
charge and statistics. 

We consider the field theory defined by the partition 
function 



2?[V',-0]exp [i (TxC) , 



(Dl) 
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We recognize Eq. (3.1) whereby contraction with the 4x4 
dimensional gamma matrices is implied by the Feynman 
slash notation and there is no TRS-breaking (Haldane) 
mass 77. 

Following Seradjeh and Franz in Ref. ^ we perform 
the family of chiral gauge transformations 



■0 =: Xc^ 



^»e75/2g+»(C"l/2)6 



(D2) 



The parameter < C < 1 implements a choice of "parti- 
tion" in the terminology of Ref. Each chiral transfor- 
mation ( |D2| ) is singular if the phase supports vortices, 
otherwise it is a pure gauge transformation. The (classi- 
cal) transformation law of C in Eq. (Dl) under the family 
of chiral transformations (D2) is 



where 



and 



c 



(D3a) 



(D3b) 



(D3c) 



Observe that, whenever 6 supports vortices and C, 7^ 1/2, 
a physical magnetic flux has appeared where there was 
none to begin with. Thus, if we demand TRS, we must 
choose ^ = 1/2 and demand that is pure gauge. 

In the spirit of Ref. ^ for general ( we define the 
family of partition functions 



V\ 



(D4) 



and compute the effective action C^^ for the gauge fields 
a^^^ and that follows from integrating the massive 
fermions and X( to lowest order in a gradient expan- 
sion. All ultra-violet divergences induced by the integra- 
tion over the fermions can be disposed of with the help 
of the Pauli-Villars regularization scheme. The effective 
action, expressed in terms of a^, a^^, and 6, that follows 
to leading order in a gradient expansion, is 



/-off 



27rm 



'5m 



Id J 
2 ^ 



2ti 



2C-1 



dj d. 



(D5) 



where Q = sgn^^i (1 — n^/m). The effective action (D5) 
fails to capture the charge of screened quasiparticles. For 
example, the conserved induced fermionic current 



27r P 



(D6) 



which is independent of the parameter C, does not repro- 
duce the induced fermionic charge ( |5.7| ) when the axial 
gauge field screens the mass vortices. It follows from 
their result that the charge bound to screened vortices 
(in which case Og^ — ^di^O = 0) is Q = (!) instead of 
g = 1/2 as found in Refs. |^, and in Sec. 0. 

Moreover, after proper dualization of the effective ac- 
tion (D5) [this dualization must include the Higgs mass, 
i.e., the first line on the right-hand side of Eq. (D5), a fact 
that was ignored in Ref. ^ it follows that the exchange 
statistics is C dependent. This is expected in view of the 
introduction of magnetic fluxes whenever C 7^ 1/2 con- 
trary to the implicit assumption made in Ref. 17 when 
choosing C 0- 
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